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ABSTRACT 


The  structure  of  a  very  strong  shock  wave  propagating 
through  a  deuterlxim-tritium  gas  mixture  and  a  pure  tritium 
gas  is  studied.  The  temperature  behind  the  shock  wave  is 
sufficiently  high  so  that  thermonuclear  reaction  probabilities 
are  large.  The  wave  structure  is  similar  to  that  of  detona¬ 
tions  in  chemically  reacting  gases,  it  is  assumed  that  the 
characteristic  times  for  collisions  and  reactions  are  such 
that  the  von  Neumann-Zeldovich  model  of  detonations  is 
applicable;  i.e.,  the  shock  can  be  treated  as  a  viscous 
gas  dynamic  shock  followed  by  a  deflagration  wave  inside 
of  which  all  the  reactions  occur.  The  physical  and  mathe¬ 
matical  assumptions  involved  in  the  analysis  of  thermo¬ 
nuclear  shock  wave  structure  are  examined.  The  reaction 
probabilities  for  deuterium  and  tritium  fusion  reactions 
are  computed  and  the  appropriate  reaction  kinetics  equations 
are  developed.  The  effect  of  energy  losses  due  to  bremsstrah- 
lung  on  the  wave  structure  is  considered  for  a  gas  that  is 
optically  thin  to  radiation  of  all  frequencies.  The  result¬ 
ing  set  of  structure  equations  are  solved  numerically  for 
several  physically  interesting  cases.  The  neutron  flux  and 
power  output  due  to  reactions  is  calculated  for  a  shock  wave 
propagating  in  a  electromagnetically  driven  shock  tube  filled 
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with  a  mixture  of  deuterium  and  tritium.  A  power  of 

3 

1  kw/cm  is  predicted  under  specified  operating  conditions. 
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CHAPTER  I 


INTRODUCTION  AND  HISTORICAL  SURVEY 


A  shock  wave  is  a  disturbance  that  propagates  with 
supersonic  speed  with  respect  to  the  undisturbed  gas  ahead 
of  it.  It  compresses  and  heats  the  gas  through  which  it 
moves.  If  the  gas  is  capable  of  undergoing  molecular, 
atomic  or  nuclear  reactions,  the  gas  composition  may  change 
across  the  shock  front, 

A  shock  wave  in  which  exothermic  reactions  take  place 
is  called  a  detonation  (see  figure  1).  A  thermonuclear 
shock  is  defined  as  a  detonation  in  which  thermonuclear 
reactions  occur.  If  the  high  temperature  created  by  the 
shock  wave  is  maintained,  the  gas  may  continue  to  react 
until  all  the  fuel  is  consumed.  Radiation  is  an  important 
energy  transfer  mechanism  at  temperatures  for  which  thermo¬ 
nuclear  reaction  rates  are  appreciable.  If  the  gas  is  opaque 
to  radiation,  a  steady  state  may  be  established  behind  the 
wave.  If  the  gas  loses  energy  through  radiation,  however, 
no  steady  state  is  established  behind  the  wave;  the  gas 
continues  to  radiate  until  its  temperature  is  the  same  as 
that  of  its  surroundings. 

The  variation  of  temperature,  pressure,  density,  speed 
and  species  concentration  with  space  inside  the  shock  wave 


FIG.  1  DETONATION  WAVE  IN  A  NON -RADIATING  GAS 
VIEWED  IN  REFERENCE  FRAME  TRAVELLING 
WITH  THE  WAVE 
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is  called  the  shock  structure.  The  flow  variables  change 
continuously  from  one  side  of  the  wave  to  the  other.  The 
interior  of  the  wave  is  where  viscous  and  thermal  conduct¬ 
ing  processes  are  important  and  where  all  the  reactions 
take  place.  However,  for  a  gas  in  which  radiation  energy 
loss  is  negligible  the  conditions  behind  the  wave  are 
independent  of  the  detailed  transport  and  reaction  processes 
occurring  in  the  shock  structure.  If  radiation  energy  loss 
cannot  be  neglected,  the  conditions  behind  the  wave  depend 
on  the  physical  processes  occurring  in  the  structure. 

In  this  paper  the  structure  of  detonations  in  gases 
capable  of  undergoing  thermonuclear  reactions  is  investi¬ 
gated.  There  is  no  known  literature  that  deals  directly 
with  thermonuclear  shock  structure.  However,  studies  of 
gas  dynamic  shocks,  detonations  in  chemically  reacting  gases, 
and  shocks  in  radiating  gases  are  relevant. 

The  early  theoretical  work  on  gas  dynamic  shock  waves 
1  2 

was  done  by  Rankine  and  Hugoniot .  They  d^r  ived  the  jump 
conditions  across  the  wave  and  considered  the  gas  dynamic 
shock  structure  equations  in  an  inviscid  gas  with  non-zero 
theritval  conductivity.  In  general,  the  gas  dynamic  shock 
structure  equations  consist  of  two  coupled  non-linear  first 
order  differential  equations.  If  either  viscosity  or  thermal 
conductivity  is  neglected,  one  of  the  equations  is  algebraic. 
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Rayleigh  studied  shock  structure  in  a  viscid  gas  with  zero 
thermal  conductivity.  The  solution  of  the  shock  structure 
equations  in  a  viscid-thermal  conducting  gas  was  obtained 

4 

by  Gilbarg  and  Paolucci  using  continuum  theory  and  by 

5  6  7  8 

Wang  Chang  ,  Mott -Smith  ,  Zoller  and  Grad  using  the 

Boltzmann  equation. 

9 

Detonations  were  first  observed  in  1881  by  Berthelot. 

Chapraan^^  and  Jouguet^^  studied  detonations  theoretically, 

postulating  that  the  gas  behind  the  detonation  moved  relative 

12  13 

to  the  wave  with  the  speed  of  sound.  Von  Neumann  ,  Dbring  , 
14 

and  Zeldovich  analyzed  the  detonation  structure  problem 
in  1940  assuming  that  a  detonation  wave  consists  of  a  gas 
dynamic  shock  followed  by  a  deflagration,  a  wave  in  which 
exothermic  reactions  take  place  that  propagates  with  sub¬ 
sonic  speed.  Under  this  assumption,  the  shock  structure 
equations  include  one  algebraic  equation  and  one  first 

order  non-linear  differential  equation.  Hirschfelder  and 
15 

Curtiss  ,  assuming  that  the  gas  dynamic  shock  and  defla¬ 
gration  are  coupled,  solved  two  simultaneous  non-linear 
differential  equations  for  the  detonation  structure.  Deto¬ 
nations  have  been  studied  for  many  kinds  of  reacting  gas 

mixtures.  For  example,  Resler  and  Cary^^  have  studied 

17 

detonations  in  dissociating  air  and  Petchek  and  Byron 
have  investigated  ionizing  shocks  in  argon. 


5 


Radiative  shock  wave  structure  has  been  studied  exten- 

18 

sively  since  1952  when  Prokof'ev  investigated  the  case  of 

steady  flow  with  zero  viscosity  and  thermal  conductivity, 

19 

i.e.,  the  "radiation  smoothed"  case.  Heaslet  and  Baldwin 

extended  earlier  work  on  the  "radiation  smoothed"  shock  by 

20 

considering  cases  with  discontinuous  profiles.  Marshak 
21 

and  Traugott  have  included  the  effects  of  viscosity  and 
thermal  conductivity  in  their  analyses  of  radiative  shock 

structure.  Shocks  in  optically  thick  atmospheres  have  been 

22  23 

investigated  by  Koch  for  a  plasma  and  Scala  and  Sampson 

for  a  chemically  reacting  gas.  Shocks  in  optically  thin 

24 

atmospheres  have  been  studied  by  Gross  for  an  ordinary 

23 

gas  and  by  Scala  and  Sampson  for  a  gas  undergoing  chemical 
reactions . 

In  the  present  paper,  the  model  of  detonation  structure 
that  von  Neumann  used  in  his  study  of  shocks  in  chemically 
reacting  gases,  is  applied  to  an  analysis  of  detonations  in 
gases  capable  of  undergoing  thermonuclear  reactions.  First, 
a  set  of  shock  structure  curves  showing  the  effect  of  reac¬ 
tions  is  obtained  for  detonations  in  tritium  and  deuterium- 
tritium  mixtures.  We  neglect  the  effects  of  viscosity  and 
thermal  conductivity  of  the  gas  and  radiative  energy  loss 
In  this  initial  treatment.  Then,  the  more  complicated 
problem  of  a  coupled  detonation  and  radiative  shock  wave 


in  tritium  and  deuterium-tritium  mixtures  is  studied.  As 
there  are  no  applied  magnetic  fields  in  the  problem, 
bremsstrahlung  is  the  only  type  of  radiation  included  in 
the  analysis.  The  gas  is  assumed  to  be  optically  thin  to 
radiation  of  all  frequencies.  A  set  of  structure  curves, 
showing  the  coupled  effects  of  reactions  and  radiation  on 
the  shock  structure,  is  obtained. 
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CHAPTER  II 

PHYSICAL  AND  MATHEMATICAL  FORMUIATION  OF  THE  PROBLEM 

In  this  chapter  the  problem  of  detonation  wave  structure 
in  gases  undergoing  thermonuclear  reactions  will  be  formu¬ 
lated  physically  and  mathematically.  Before  writing  down 
the  appropriate  shock  structure  equations  for  thermonuclear 
shocks  (first,  assuming  no  radiation  losses,  and  then  includ¬ 
ing  radiation  losses),  the  von  Nevimann  model  of  detonation 
structure  will  be  described.  A  discussion  of  the  physical 
and  mathematical  assumptions  used  to  simplify  the  structure 
equations  will  follow. 

2.1  Von  Neumann  Model  of  Detonation  Structure 

The  model  of  detonation  structure  used  by  von  Neumann 
in  his  study  of  shocks  in  chemically  reacting  gas  mixtures, 
treats  a  detonation  as  a  shock  followed  by  a  deflagration. 

The  shock  wave,  which  is  thin  compared  to  the  deflagration 
that  follows  it,  propagates  supersonically  into  the  undis¬ 
turbed  and  unreacted  gas ;  the  deflagration,  inside  of  which 
all  the  reactions  occur,  propagates  subsonically  into  the 
gas  through  which  the  shock  has  already  travelled. 

The  changes  in  pressure,  density,  speed,  temperature 
and  ratio  of  reaction  product  density  to  total  density  (  ) 

inside  the  detonation,  viewed  in  a  reference  system  moving 


with  the  detonation^  are  depicted  schematically  in  figure  2. 
Eadi  variable  asymptotically  approaches  its  upstream  value  as  x 
goes  to  -  GP  and  its  downstreaun  value  as  x  goes  to  +  oo  , 

In  the  shock  wave  (points  1  to  2  in  figure  2),  the  gas  is 
ccxnpressed,  decelerated  and  heated.  In  the  deflagration 
(points  2  to  3  in  figure  2),  the  gas,  as  the  reactions 
proceed,  is  gradually  expanded  and  accelerated.  The  tempera¬ 
ture,  which  reaches  almost  its  final  value  behind  the  shock, 
rises  to  a  maximum  inside  the  deflagration  and  then  decreases 
slightly  as  the  fuel  gets  used  up  and  the  random  thermal 
energy  is  transferred  into  bulk  kinetic  energy  of  the  gas. 

The  von  Neumann  model  is  a  valid  description  of  a 
detonation  provided  thal:  the  characteristic  time  for  reac¬ 
tions  to  take  place,  is  long  compared  to  the  time  involved 

25 

in  passage  of  the  shock  wave  (Hirschfelder,  p.  801).  The 

time  for  the  shock  wave  to  pass  in  a  plasma  depends  on  the 

26 

ion-ion  collision  time,  behind  the  shock  (Von  Karman  ) . 
Therefore  the  validity  criterion  for  the  von  Neumann  model 
is  that  <  <  \  behind  the  gas  dynamic  shock. 

In  order  to  get  some  physical  feeling  for  the  conditions 
inside  thermonuclear  shocks,  we  shall  note  some  typical  values 
of  parameters  such  as  temperature,  pressure  and  speed.  The 
initial  gas  mixture  (point  1  on  figure  2)  is  assumed  to  oe 
fully  ionized  to  avoid  the  mathematical  complications  of 
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....ocution  .nd  ionization.  «.aro,an  ia  fully  ionized  at 

rt^utSzlO^OK  (Btezing”).  Initial  nu.*er  deneitie. 

of  intereat  for  controlled  «.er«=nuclear  fuaion  are 

^  io«  -  io“  ^ 

In  apeed.  WHiC  ia  alao  t.e  initial  ga.  apeed  relative  to 
a  reference  system  moving  with  the  detonation,  is  abo 

lO^cm/sec. 

«,e  ahoc.  vave  heata  the  gaa  to  a  ten^rature  of  about 
OK.  Behind  it  Ipoint  2  on  figure  2),  nu^r  denartrea 
,.ve  increased  by  a  factor  of  four  and  gas  speeds  have 
decreased  by  the  sane  factor  (see  section  2.3).  Aa  the 
reactions  go  to  completion  inside  the  deflagration  wave, 
the  density,  speed,  pressure  and  temperature  change  no  mo 

than  a  factor  Of  two.  ihe  characteristic  temperature,  and 

nomber  densities  behind  the  deflagration  (point  3  on  igure 

sre  about  10^“  °K  and  10^=  o.'^  respectively. 

«,e  temperature,  and  speed,  characteristic  of  thermo- 

„oclear  ahoch.  are  much  higher  than  those  a.aociated  with 

other  hind,  of  strong  ahocK  waves.  A  graph  of  Poat-ahoCh 

■  -.  for  hydrogen  is  shown  in  fiBure  3. 
temperature  vs.  wave  speed  for  hyorog 

,-sd  is  increased,  the  gas  first  becomes  disso- 
As  the  wave  speed  is  mere  »  , 

Wed  As  the  wave  speed  is  increased 
ciated  and  then  ionized.  As  me 

sioes  in  the  gas  become  iimportant. 
f^^ther,  radiation  processes  in  the  g 

•«  rtiven  by  the  Fankino-Hugoniot 
The  post-shock  temperature  is  given  by  tn 
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conditions  for  a  gas  dynamic  shock.  In  the  dis^sociation 
and  ionization  speed-temperature  regimes,  the  appropriate 
chemistry  has  been  incorporated  into  the  jiunp  conditions 
(Taussig  ) .  In  optically  thick  atmospheres,  the  radiation 
energy  density  becomes  important  compared  to  the  internal 
energy  of  the  gas  for  >  lO^cm/sec  (Koch^^) .  Radiation 
terms  must  be  included  in  the  jump  conditions  for  such  shock 
speeds.  For  an  optically  thin  atmosphere,  radiation  energy 
losses  become  important  only  for  shock  speeds  above  lO^cny^sec. 
There  are  no  Rankine-Hugoniot  conditions  for  shocks  in  opti¬ 
cally  thin  atmospheres  because  no  steady  state  is  established 
behind  the  shock.  The  post-shock  temperature  shovm  in 
figure  3  is  the  maximum  temperature  reached  in  the  structure 
of  the  shock  (Gross^^) . 

The  post-shock  temperature  for  thermonuclear  shocks  is 
the  temperature  behind  the  entire  wave,  gas  dynamic  shock 
plus  deflagration.  It  is  given  by  the  set  of  Rankine- 
Hugoniot  conditions  that  include  the  energy  released  and 
the  change  in  gas  composition  due  to  reactions  (see  section 
2.3).  The  lowest  detonation  speed  which  satisfies  these 
Rankine-Hugoniot  conditions  is  called  the  Chapman-Jouguet 
speed.  If  the  detonation  travels  with  the  Chapman-Jouguet 
speed,  the  downstream  flow  is  sonic.  For  speeds  greater 
than  the  C-J  speed,  there  are  two  solutions  to  the  Rankine- 


13 


Hugoniot  conditions.  The  upper  branch,  a  strong  detonation, 

has  subsonic  flow  behind  it;  the  lower  branch,  a  weak  deto- 

30 

nation,  has  supersonic  flow  behind  it  (Hayes  ) .  A  graph 

of  post-shock  temperature  as  a  function  of  wave  speed  for 

deuterium  is  shown  in  the  figure  3.  Deuterium  and  tritium 

fusion  reactions  have  much  larger  cross  sections  than 

hydrogen  fusion  reactions  and  are  therefore  of  greater 

interest  in  problems  of  controlled  thermonuclear  fusion 

31 

(Taylor  and  Tobolsky  ) . 

2.2  Physical  and  Mathematical  Assumptions 

The  physics  to  be  described  in  a  complete  treatment  of 
a  thermonuclear  shock  structure  problem  includes  (1)  reac¬ 
tions,  v2)  radiation,  (3)  transport  processes,  (4)  non-equi¬ 
librium  effects  due  to  the  multi-component  nature  of  gas, 
and  (5)  the  effects  of  applied  electromagnetic  fields.  Reac¬ 
tions  occur  in  the  gas  which  has  been  shock  heated.  Some 
of  the  energy  added  to  the  gas  by  reactions  is  lost  through 
radiation.  Transport  processes  due  to  the  viscosity  and 
thermal  conductivity  of  the  gas  account  for  the  bulk  flow 
of  momentum  and  energy.  Since  the  gas  is  a  mixture  of 
several  kinds  of  ions  and  electrons,  non-equilibrium  pro¬ 
cesses  (e.g.,  diffusion,  charge  separation,  and  energy  and 
momentum  transfer  between  species)  occur.  External  magnetic 


fields,  necessary  to  contain  a  thermonuclear  plasma  in  order 
to  tap  its  energy,  further  complicate  the  physics  by  adding 
cyclotron  radiation  and  hydromagnetic  effects. 

A  mathematical  treatment  of  thermonuclear  shock 
structure  including  all  the  relevant  physics  would  require 
the  solution  of  the  set  of  non-linear  differential  equations 
that  result  from  a  coupling  of  the  many  species  Boltzmann 
or  continuum  equations  with  the  radiative  transfer  equation 
and  Maxwell's  equations.  To  reduce  the  number  of  differen¬ 
tial  equations  in  the  problem  to  a  manageable  set,  some 
of  the  physical  effects  must  be  igr.ored.  The  present  treat¬ 
ment  will  include  reactions  and  radiation  while  neglecting 
transport  phenomena,  non-equilibrium  effects  due  to  the 
presence  of  more  than  one  species,  and  the  containment  problem. 
The  many  species  nature  of  the  gas  is  incorporated  in  the 
state  and  continuity  equations  only. 

A  list  of  assumptions  used  in  the  treatment  of  thermo¬ 
nuclear  shock  structure  without  radiation  in  tritixim  and 
deuterium-tritium  mixtures  is  given  below.  It  is  followed 
by  a  discussion  of  their  j^ysical  validity.  The  further 
assumptions  that  must  be  made  when  radiation  is  included 
in  the  problem  are  then  listed  and  discussed. 

PHYSICAL  ASSUMPTIONS  FOR  THERMONUCLEAR  SHOCK  STRUCTURE 


1.  Ihe  von  Neumann  model  of  detonations  is  valid: 


vf' 


15 


the  shock  wave  is  thin  compared  to  the  deflagration 
wave;  all  reactions  take  place  in  the  deflagration 
region.  Since  the  influence  of  reactions  on  the 
detonation  structure  is  of  primary  interest,  the 
shock  is  treated  as  a  discontinuity  and  only  the 
deflagration  structure  is  studied. 

2.  Transport  phenomena  due  to  viscosity  and  thermal 
conductivity  are  negligible  in  the  deflagration. 

3.  Non -equilibrium  effects  due  to  the  presence  of 
more  than  one  species  are  negligible:  All  species 
travel  with  the  local  gas  speed;  the  particles  of 
each  species  have  a  Maxwellian  distribution  of 
energies  at  the  local  gas  temperature;  charge 
separation  is  negligible. 

4.  One-dimensional,  single  fluid,  continuum  equations 
provide  an  adequate  description  of  the  phenomena. 

5.  All  macroscopic  velocities  and  temperatures  are 
low  enough  so  that  relativistic  effects  can  be 
neglected. 

6.  There  are  no  applied  electromagnetic  fields. 

Assumptions  1,  2  and  3  were  made  in  order  to  simplify 

both  the  mathematics  and  the  physics  of  the  thermonuclear 
shock  structure  problem.  Their  physical  validity  depends 


on  the  ratios  of  characteristic  collision  and  reaction  times 
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for  typical  temperatures  and  number  densities  inside 

thermonuclear  shocks. 
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Spitzer  defines  the  following  characteristic  times 
for  a  plasma:  The  time  required  for  a  single  species,  non- 
Maxwellian  distribution  to  approach  the  Maxwellian  distri¬ 
bution  is 

sec,  (2.1 

n(cm“^)  jin  A 


where  Z  is  the  atomic  number  of  tlie  ion,  A  is  its  atomic 
weight,  and 


./V  = 


k^  ^ 


2  e^ 


IT  n  / 
e  / 


(2.2) 


The  self  collision  time  for  electrons,  X  ,  is  just  'C 

ce  c 

with  A^=  1/1836,  the  atomic  weight  of  electrons.  The  rate 
at  which  equipartition  of  energy  between  two  groups  of 
particles  is  established  is 


X 


5.87 


Ai  A^ 


2  2 

n^Zj^Z^  xnA 


(2.3) 


If  the  two  groups  are  ions  and  electrons, 
then 


X  = 

eq 


A 

5 .87  X 
A* 


T  *  (°K) 

n  z2  in  A 


e 


and  if 


■'2' 


sec. 


(2.4) 
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The  time  associated  with  viscous  processes,  which  is  what 
26 

von  Karman  calls  the  characteristic  time  associated  with 
thermodynamic  transformations,  is  the  ion -ion  collision 
time.  The  time  associated  with  thermal  conduction  processes 
is  the  ion-electron  collision  time.  Another  relevant 
characteristic  time  is  the  mean  time  between  reactions 


T  =  _i _ 

R  n  <c  V  >  * 


(2.5) 


where  <.rv  >  is  the  reaction  cross  section  times  the  relative 

speed  of  the  interacting  particles  averaged  over  the  product 

35 

of  their  distribution  functions  (Osborn  ) . 

The  ratios  of  the  collision  times,  X  ,  t  ,  L  are 

*  c*  ce*  eq 

independent  of  temperature  and  number  density.  For  a  tritium 
plasma  (A  ■  3) , 

^  =  74.3  and  =  38.2. 

^ce  Tc 

Ratios  of  collision  times  to  the  reaction  time  must  be  eval¬ 
uated  separately  for  each  temperature  and  number  density 
of  interest.  For  n  =  10^^  cm  ^  and  10^  ®K  <  T  <  10^®  ^K, 

In  A.  (which  is  not  a  very  sensitive  function  of  temperature 

and  number  density)  ^21.  The  values  of  X  eq/  Cr  and 
T  T 

^c/  R  (with  A  =  3)  are  given  in  Table  1  for  several 
temperatures  in  this  range. 
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T  (Ok) 

Teq/tR 

Tc/Tr 

9 

1  X  10 

-17 

2  X  10 

0.02 

0.0006 

9 

5  X  10 

-17 

8  X  10 

1.02 

0.027 

1  X  10^° 

,  i 

1  X  10  i 

I 

3.8 

0.11 

2  X  10^° 

3  X  10“^®  ; 

t 

30.5  j 

t 

0.79 

TABLE  1.  Ratios  of  Characteristic  Times 


The  shock  region  is  thin  compared  to  the  deflagration 

region  (assumption  1)  if  Tc/Tr  <.cl  behind  the  shock. 

For  a  tritium  plasma  the  temperature  inside  the  shock  rises 

4  10 

from  its  downstream  value,  5  x  10  °K,  to  10  before  any 

reactions  take  place.  However,  as  the  temperature  approaches 
its  post-shock  value,  ion-ion  collision  times  and  reaction 
times  become  comparable.  Therefore  the  reactions  may  not 
be  completely  confined  to  the  deflagration  wave;  some  may 
take  place  in  the  tail  of  the  shock. 

Viscosity  is  negligible  inside  the  deflagration  (assump¬ 
tion  2)  if  I^c/Cr  1  for  typical  temperatures  and  number 
densities.  The  requirement  for  neglecting  thermal  conduc¬ 
tion  and  diffusion  is  that  l/H^  ^  ^  Karman^^) 
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M  is  the  Mach  number,  the  ratio  of  the  local  gas  speed  to 


the  local  sound  speed: 


M  = 


2  .i 


iR  T  /  * 


(2.6) 


y  ,  the  ratio  of  specific  heats  at  constant  pressure  to 
that  at  constant  volvime,  has  the  value  5/3  for  a  plasma, 
a  mixture  of  monatomic  gases.  It  can  be  shovm  (see  section 
2.3)  that  inside  the  deflagration. 


Since 


^  ^  4  1. 


c/  L  R  )  (^  l/M^j  > 


1  inside  the  deflagration,  the 


effects  of  viscosity,  thermal  conduction  and  diffusion  would 
have  to  be  included  in  a  fuller  treatment. 

The  requirement  for  the  ions  and  electrons  to  be  in  ther¬ 
mal  equilibrium  (assumption  3)  is  that  ^eq/'I^R  c  1  inside 
the  deflagration.  The  reaction  product  ions  are  released  with 
Mev  energies.  The  electrons,  since  they  are  not  involved  in 
reactions,  are  cooler  than  the  ions.  Although  reaction 
product  ions  transfer  energy  to  fuel  ions  in  a  time  comparable 
to  the  reaction  time,  energy  transfer  from  the  hot  ions  to 
the  cooler  electrons  occurs  much  more  slowly.  (For 
T  =  2  X  10^^  OK,  C  c/  t  R  '  0.8  while  'Xeq/^R  30)  . 

The  neutrons  released  in  reactions,  (which  we  assume  are 


^  - 


'•  ■  ■  A' 


.  I 


llBBSiiiMn^iiiti » i  . — - Jl_,_ 


■ikb.^  t  wr-  fe  •  T>.  V-SW  ■!•.•'  • 


20 


kept  in  the  gas),  must  be  in  equilibrium  with  the  ions  and 

electrons  for  the  gas  to  be  in  thermal  equilibrium.  Cross 

sections  for  neutron-triton  scattering,  t:,,  are  w  1  barn 

N 

33 

for  Mev  neutron  energies  (Hughes  and  Schwartz  ) .  A  10 

Mev  neutron  goes  through  2-5  collisions  before  its  energy 

34 

reaches  2  Mev  (Glasstone  ) .  The  neutron  collision  time, 

Xc«  (number  of  collisions  to  thejrmalize  )  x  l/(nir,v„). 

"  N  N 

10  —  9 

For  T“10  OK,  V„  ia  10  cro/sec.  Therefore  the  ratio  of 

N 

the  neutron-ion  collision  time  to  the  ion-ion  time,trcjj/'Cc  is^l, 
and  the  neutrons  therroalize  with  the  ions  in  a  time 
comparable  to  the  ion-ion  collision  time.  Equilibration 
times  for  neutrons  and  electrons,  however,  are  much  larger 
than  those  for  neutrons  and  ions  since  neutron-electron 
collision  cross  sections  are  very  small.  The  neutrons, 
like  the  ions,  are  therefore  not  in  equilibrium  with  the 
electrons . 

The  gas  has  local  charge  neutrality  if  the  electron 
Plasma  frequency,  , 


r 


4  -n  n  e 
e 

m 


2  \ 


(2.7) 
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(Spitzer  )  is  much  larger  than  the  frequency  associated 
with  the  shocK  wave,  typical 

temperatures  and  number  densities  mentioned  previously. 


1 
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Iv 

R 

Therefore  the  assumption  of  charge  neutrality  is  a  good  one. 

One -dimensional,  single  fluid,  continuum  equations 

(assumption  4)  have  been  used  in  many  studies  of  shock 

4  22 

structure  (e  .g .  Gilbarg  and  Paolucci  and  Koch  ) .  A 

discussion  of  the  use  of  continuum  fluid  equations  rather 

than  the  Boltzmann  equation  is  given  by  Gilbarg  and 
4  24 

paolucci  and  by  Gross,  Treatment  of  the  ions,  neutrons 
and  electrons  as  a  single  fluid  neglects  the  lack  of  thermal 
equilibrium  between  the  electrons  and  the  ions  and  neutrons. 

If  the  electrons  are  assumed  to  be  as  hot  as  the  ions  and 
neutrons,  a  significant  fraction  of  them  would  have  relati¬ 
vistic  energies.  The  non-relativistic  equations  are  valid 
(assumption  5)  if  kT  <<  me  for  each  species,  (me  for 
electrons  is  ~  |  Mev.)  The  use  of  non-relativistic  equa¬ 
tions  is  justified  physically,  however,  since  the  electrons 
are  actually  cooler  than  the  ions  and  neutrons. 

The  following  physical  assumptions  are  made  in  addition 
to  those  listed  previously  when  radiation  is  included  in 
the  thermonuclear  shock  structure  problem: 

7.  Radiation  effects  are  significant  in  the  deflagration 


region  only. 


’ .  -  ■  ,V 
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8.  The  gas  is  optically  thin  to  radiation  of  all 
frequencies . 

9.  The  radiation  comes  from  ion-electron  bremsstrah- 
lung  only. 

Since  the  von  Neiunann  model  of  detonations  is  assumed, 
the  effects  of  radiation  are  studied  in  the  deflagration 
region  only.  In  figure  3,  it  is  seen  that  for  thermonuclear 
shock  speeds,  the  maximum  temperature  behind  an  optically 
thin  shock  is  only  slightly  lower  than  the  temperature 
behind  a  non-radiative  gas  dynamic  shock.  Therefore  assump¬ 
tion  7  is  valid. 

A  gas  is  optically  thin  to  radiation  if  the  photon 

absorption  mean  free  path,  frequencies,  V  , 

is  much  larger  than  all  the  other  cdiaracteristic  lengths  in 

the  problem.  The  reaction  mean  free  path  is  the  largest 

characteristic  length  for  non-radiative  thermonuclear  shocks. 

Therefore  the  gas  is  optically  thin  if  \  >  >  X_  for 

P 

all  frequencies.  (For  T  =  10^^  OR  and  n  =  10^^  an 

9 

is  ~  10  an) .  The  only  type  of  radiation  occurring  in  a 
fully  ionized  plasma  in  the  absence  of  magnetic  fields  is 
bremsstrahlung .  Since  non-relativistic  equations  are  used, 
only  ion-electron  bremsstrahlung  is  considered;  electron- 
electron  bremsstrahlung  is  neglected.  An  approximate  expression 
for  the  mean  free  path  for  absorption  of  a  photon  by  a  hydrogen 


f 


'•V 


isotope  plasma,  as  the  result  of  a  free-free  transition, 


is  (Glasstone  and  Lovberg  ), 


^  ii«“5  (hev)  ; 

7  X  10  - s - -  cm. 


(2.8) 


The  power  emitted  per  unit  wavelength  per  unit  volume  due 
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to  ion-electron  bremsstrahlung  (Glasstone  and  Lovberg  ) 


11^6.01  X  10“^\  ^2)  T-^X"^exp(-^^)^£|5~.  (2.9) 

^  ”  cm'^  A-sec 


The  maximum  value  of  occurs  for 

dA 


-  8. 

T  (kev) 


(2.10) 


For  T  =  1  Mev,  X  is  6.2  x  10  ^  X  (V  ^  5  x  10^®  cps). 
'  max  max 


The  photon  absorption  mean  free  path  for  this  temperature 
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and  frequency  is  ~  10  cm,  which  is  much  larger  than  the 


reaction  mean  free  path.  The  smallest  frequency  for  which 


is  >  >  X  p  is  -  lO"^  cps  (  ^  »  3  X  lO"^  A) .  The 


ratio  of  the  power  emitted  at  a  frequency  of  10  cps  to 


20  -11 

that  emitted  at  a  frequency  of  5  x  10  cps  is  3  x  10 


Therefore  the  fraction  of  radiation  energy  emitted  at 


frequencies  for  which  the  gas  is  not  optically  thin,  is 


negligible.  It  is,  therefore,  a  good  assumption  to  treat 
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the  gas  as  an  optically  thin  atmosphere. 

2.3  Shock  Structure  Ecmations  Without  Radiation 

The  following  are  the  time-*dependent  continuum  fluid 
conservation  equations  for  the  deflagration  wave  (von 
Karman^^) : 

till 

The  continuity  equation  for  the  i^  specie  is 

If  i  (Hi  u)  -  W. ,  i  =  1,2,...,  N,  (2.11) 

where  N  is  the  nuitiber  of  species  in  the  mixture,  and  is 

th 

the  production  rate  of  the  i  species (see  section  2.4). 
Multiplying  each  continuity  equation  by  m^  and  summing  over 
all  species,  one  obtains  the  mass  flow  equation, 

^  ^  (f^u)  =  m.  W  ,  (2.12) 

where  n^  m^  .  (2.13) 

If  no  particles  (e.g.,  neutrons)  escape  from  the  gas,  and 
if  negligible  mass  is  lost  or  created  in  reactions,  then 
m.  W.  «  0.  The  mass,  ^m,  lost  in  a  nuclear  reaction  is 

‘Y  1  i  *  * 

/^im  «  (mj^  +  m^)  -  (m^  +  m^)  ,  (2.14) 

where  m^  and  m^  are  the  masses  of  the  reacting  nuclei  and 
m^  and  m^  the  masses  of  the  product  nuclei  (Evans^®) .  If 


Am/(ra2^+»n2)  is  <  <  1,  the  mass  lost  in  reactions  can  be 
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neglected  in  the  mass  equation.  For  thermonuclear  reac¬ 
tions,  +  m^)  is  <^0.01,  so  mass  is  conserved  to  a 

high  degree  of  accuracy. 

A  momentum  flow  equation  for  each  species  in  the  mixture, 
analogous  to  the  species  continuity  equations  for  mass  flow, 
can  also  be  written  down.  The  momentum  flow  equation  for 
the  mixture  is  obtained  by  summation  of  the  species  momentum 
equations.  If  the  viscosity  and  thermal  conductivity  of 
the  mixture  and  the  diffusion  of  the  species  are  neglected, 
the  momentum  flow  equation  for  the  mixture  is, 

^  {^u)  (pu^)  =.|£  ,  (2.15) 

where  p  =  £  k  T  .  (2.16) 

The  energy  flow  equation  for  the  mixture,  which  is  derived 
by  summing  the  energy  e(^ations  for  all  the  species,  is 

-ft  -  ‘’ra»  ■ 

where  e  is  the  internal  energy  of  the  mixture,  and  is 
a  volume  energy  source  due  to  thermonuclear  reactions. 

If  the  shock  is  viewed  frcmi  a  coordinate  system  moving 
with  the  shock  speed,  all  time  derivatives  in  the  continuity, 
mass,  momentum  and  energy  equations  vanish.  If  the  result¬ 
ing  time-independent  equations  are  integrated  once  with 
respect  to  the  spatial  variable,  x,  they  become 
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pn  =  A,  (2.18) 

p  +  pu^  =  B  ,  (2.19) 

2 

and  ^  1  '  ^ 

where  h  =  pu  +  e  =  ^  T  is  the  enthalpy  of  the  gas 

(which  is  assumed  to  be  ideal) . 

The  constants  of  integration.  A,  B,  and  C, evaluated 
from  the  steady  state  conditions  upstream  and  downstream 
from  the  shock  (see  figure  2),  are 


u. 


^  ^3  ^3  * 


(2.21) 


B  =  -  Pj  +  /3  , 

2  2 

and  C=h^+^  =h3+l!3  -  Q  • 

2  2 


(2.22) 

(2.23) 


In  state  1  the  reactions  have  not  begun  so  =  0.  The 
reacx-.  ons  have  gone  to  completion  in  state  3  so  =  Q, 
the  energy  released  per  unit  mass  of  the  initial  mixture 
if  it  reacts  completely. 

Thermonuclear  shocks  are  very  strong  shocks,  i.e., 

2 

Pj^  '  ^  Rankine-Hugoniot  conditions  for  a  thermo 

nuclear  shock  follow  from  equations  (2.21),  (2.22),  (2.23) 
and  the  assumption  of  a  strong  shock.  The  ratio  of  the 
gas  speed  behind  the  thermonuclear  shock  to  the  gas  speed 
ahead  of  the  thermonuclear  shock,  (the  shock  speed) ,  is 
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J_  (l-  0/u,^)" 

“l  ‘  ft  ‘  1 

where  ^  =  |  . 

The  minimum  shock  speed  that  satisfies  equation 
Chapman -Jouguet  speed,  is 

j. 

u^  =  C2(  jr^  -  1)  of. 


(2.24) 


(2.24),  the 


(2. 25) 


For  C-J  flow  conditions,  equation  (2.24)  reduces  to 


y  5 

u^  “  F+r  “  8 


(2.26) 


The  ratio  of  the  gas  speeds  behind  the  gas  dynamic  shock  to  that 
in  front  of  it,  ^  ,  is  given  by  equation  (2.24)  with  Q  =  0. 


The  two  values  of  ^  that  satisfy  equation  (2.24)  with 

^1 

Q  =  0  are 


jr-1 

1 

1 

and  ^ 

^1 

^  r+T 

4 

u 

(2.27) 


One  feature  of  thermonuclear  shocks  is  that  the  tempera¬ 
ture  of  the  gas  behind  the  thermonuclear  shock,  T^^  is  indepen¬ 
dent  of  the  initial  temperature  of  the  gas.  T^  is  related 
to  the  energy  released  in  reactions,  Q,  and  the  shock  speed. 


by 
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1-2  (IT  -1)  fiL 


'1  -1 


u. 


(2.28) 


For  C-J  flow. 


R-T 

W3 


(m  ) 


y  Q 


5 


Q. 


(2.29) 


The  temperature  of  the  gas  behind  the  gas  dyneunic  shock, 
T^,  is  given  by  equation  (2.28)  with  Q  *=  0,  i.e.. 


R  T,  2(if-l)  2 

—  «  -  ■  -  u  • 

W2  (ir+l)2  1 


If  the  shock  is  propagating  with  the  C-J  speed, 

RT  (iT-l)^  2 

=  4  -  Q  =  f  Q  .  (2.31) 

^2  Y+l  ^ 


The  temperature  behind  the  gas  dynamic  shock,  T^,  and 

that  behind  the  thermonuclear  shock,  T^,  are  plotted  as 

functions  of  the  shock  speed  in  figure  4.  The  gas  is  assumed 

to  be  a  deuterium-tritivim  mixture  containing  equal  numbers 

of  tritixim  and  deuterium  ions.  The  C-J  speed  for  the  mixture 
q 

is  3.48  X  10  cm/sec  corresponding  to  a  post-detonation 
temperature  of  4.2  x  10^®  °K.  T^  is  higher  them  T^  for  all 
shock  speeds.  The  increase  in  specific  enthalpy  across  the 
deflagration  is  accompanied  by  an  increase  in  temperature, 
since  the  mean  molecular  weight  of  the  mixture  remains  con¬ 
stant  as  the  reactions  occur. 


FIG.  4  POST-SHOCK  AND  POST  DETONATION  TEMPERATURE 
VS.  DETONATION  SPEED  FOR  DEUTERIUM-TRITIUM  MIXTURE 

(hi) 


FIG.  5  POST-SHOCK  AND  POST-DETONATION  TEMPERATURE 
VS.  DETONATION  SPEED  FOR  TRITIUM 
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The  temperatures  behind  the  gas  dynamic  and  thernra- 
nuclear  shocks  as  a  function  of  shock  speed  for  a  tritium 

9 

gas  are  shovm  in  figure  5.  The  C-J  speed  is  2.5  5  x  10  ci^/sec, 
creating  a  temperature  of  2,2  x  10^®  OK  behind  the  thermo¬ 
nuclear  shock.  The  temperature  behind  the  deflagration 
(thermonuclear  shock)  is  equal  to  the  temperature  behind 
the  gas  dynamic  shock  for  shocks  propagating  with  the  C-J 
speed,  T^  is  T^  for  shocks  propagating  with  speeds 

9 

between  the  C-J  speed  and  3,4  x  10  ciV'sec,  If  the  shock 
speed  is  increased  further,  the  temperature  behind  the 
deflagration  will  be  lower  than  that  behind  the  gas  dynamic 
shock.  An  increase  in  the  enthalpy  of  the  gas  as  the  reac¬ 
tions  occur  does  not  necessarily  show  up  in  increased 
temperatures,  as  the  mean  molecular  weight  of  the  mixture 
decreases  with  reactions.  The  temperature  behind  a  gas 
dynamic  shock  propagating  in  a  gas  containing  the  reaction 
products  of  tritium  fusion  reactions,  (alpha  particles  and 
neutrons  in  a  1:2  ratio  by  number),  is  also  shown  in  figure 
4  for  comparison. 

2,4  Reaction  Rates 

The  source  terms,  W^,  in  the  species  continuity  equations, 
must  be  specified  in  order  to  solve  the  non-radiative  shock 
structure  eqvaations  across  the  deflagration.  These  functions 
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are  related  by  the  requirement  for  mass  coiisemration, 

^  mi  Wf  =  0.  It  can  be  shovm  that  they  all  depend  on 
the  reaction  rates  for  the  fusion  reactions  taking  place 
in  the  plasma. 

The  reaction  rate  for  binary  nuclear  reactions  is 


w 


”l  ”2 


1+  X 


(<rv> 


(2.32) 


where  n^^  and  n^  are  the  number  densities  of  the  interacting 
nuclei  and  where  the  Kronecker-  S  takes  into  account  the 
possibility  that  the  reacting  particles  are  of  the  same 
kind  =  A^) .  If  each  species  has  a  Maxwellian  distri¬ 
bution  of  velocities  at  the  same  temperature,  the  reaction 
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probability  <t  nj  >  is  given  by  (Glasstone  and  Lovberg  ) 


2  - 

V  dv, 


(2.33) 


where  yu.  =  — is  the  reduced  mass  of  the  system, 

12 

T"  =  cross  section  for  binary  nuclear  reactions,  and 

v  =  relative  speed  of  the  interacting  nuclei. 

It  is  convenient  to  express  <'iV>  as  an  integral  over 

center-of-roass  energy,  E  .  Since 

cm 

=cm  =  2/^^^  » 


(2.34) 
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the  reaction  probability  can  be  written, 


<T.>  =  in.  f_i_f  1 


■4  (E^ _ )  exp  (  -  _an 

]<T , 


E  dE  .  (2.35) 
cm  cm 


Nuclear  reactions  are  studied  in  the  laboratory  by  sending 
beams  of  nuclear  projectiles  toward  a  target.  Therefore 
nuclear  reaction  cross  sections  are  usually  given  as  a 
function  of  the  laboratory  energy  of  the  incident  particle 
E  ,  If  the  target  nuclei  are  in  motion,  as  in  a  plasma, 

E_  is  the  relative  energy  of  the  interacting  nuclei.  With 

Jj 

the  relation  between  E,  and  E  , 

L  cm 


(2.36) 


the  reaction  probability,  written  in  terms  of  E_ ,  is 

Li 


}  r 

=  §IL  / 

•P  \  . 

2 

ro  ' 

k2  n  kr 

) 

c 

reaction 

mean 

free  path. 

(  . 

j  'i'(E^)exp/- 


It  1  E  <3  E  .  (2.37) 
m  KT/  ^  ^ 


^  ,  and  a  reaction  power  density  can  be  defined  in  terms 

X\ 

of  <^'J>  .  If  we  let 


n  <'*v/>  * 


the  reaction  mean  free  path  is  defined  as 


'C 


(2.38) 


(2.39) 


1 


/p  kT\^ 

Where  ^  Maxwellian  velocity  distribution,  (2.40) 

The  power  density  is  just  the  reaction  rate  multiplied  by 
the  energy  released  per  reaction. 


‘TN 


(2.41) 


where  Q  =  (Am)  c  ' . 
R 


2.5  Shock  Structure  Equations  With  Radiation 

No  steady  state  is  established  behind  thermonuclear 
shocks  propagating  in  a  gas  that  is  losing  energy  through 
radiation.  Therefore,  there  are  no  Rankine-Hugoniot  condi¬ 
tions  which  determine  the  post-shock  state  of  the  gas  given 
its  initial  state.  Since  no  steady  state  is  established 
behind  the  shock,  it  is  not  obvious  for  the  radiative  case 
that  a  reference  frame  travelling  with  the  shock  speed  is 
one  in  which  all  time  derivatives  in  the  structure  equations 
vanish. 

A  physical  argument:  can  be  used  to  show  that  the  shock 
frame  is  a  steady  reference  frame  for  shocks  in  optically 
thin  atmospheres.  It  is  apparent  that  successive  particles 


of  gas,  passing  through  the  shock,  go  through  the  same  history 
in  any  region  containing  the  shock  that  is  much  thinner  than 
a  photon  absorption  mean  free  path.  The  profiles  of  the 
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flow  variables,  viewed  in  the  shock  frame,  are  therefore 
stea^  throughout  such  a  region.  Thus,  we  can  validly 
set  time  derivatives  equal  to  zero  in  the  structure  equa¬ 
tions  if  we  limit  the  integration  of  these  equations  to 
a  region  containing  the  shock  that  is  much  thinner  than 
a  photon  adasorption  mean  free  path.  (The  structure  equa¬ 
tions  may  be  integrated  over  all  space  for  thermonuclear 
shocks  in  a  non-radiating  gas.) 

All  of  the  non-radiative  thermonuclear  shock  struc¬ 
ture  equations  are  valid  for  thermonuclear  shocks  in 
optically  thin  atmospheres  except  the  energy  equation. 

The  radiation  pressure  in  the  radiative  momentum  equation 
and  the  energy  density  term  in  the  energy  equation  are 
negligible  because  all  of  the  radiation  is  assumed  to 

escape  from  the  gas.  The  only  radiative  term  in  the  energy 
dq 

equation  is  ^  ,  the  divergence  of  the  radiative  flux, 
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(Vincenti  and  Kruger  ) .  The  energy  flow  equation  with 
radiation  for  thermonuclear  shocks  in  an  optically  thin 
atmosphere  is 

(  pu  (h  +  ^  ^  =0.  (2.42) 

If  equation  (2.42)  is  integrated  over  all  space  and  then 
divided  by  the  mass  flow  constant,  ^  u  =  A,  then 


h  +  ~  - 


2  "  S?N  ’*' 


i  f  <iq 
A  ^  7“, 

o  dx' 


<Sx*  =  C. 


(2.43) 


An  e:q)ression  for  ^  can  be  obtained  from  the  radia¬ 
tive  transfer  equation. 


where  is  the  radiative  intensity ,ya  is  the  cosine  of 
the  angle  between  the  x-axis  and  the  direction  of  observa¬ 
tion  of  I^,  is  the  volumetric  emission  coefficient,  and 
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f-  IS  the  volumetric  absorption  coefficient  (Goulard  ) . 

For  an  optically  thin  gas,  emission  processes  are  more 

important  than  absorption  processes  so 
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(Vincenti  and  Kruger  ) .  The  radiative  transfer  equation 


thus  becomes. 


yti  -  =n.  •  (2.45) 

dx 

If  it  is  integrated  over  all  frequency  and  solid  angle, 

assuming  that  the  emission  is  isotropic,  we  get  an  expression 

for  ^  , 
dx  * 


^  =  4  j  ,  dv  , 


where  q  =  2  v  |  j  kI  d  w  dy  . 


(2.46) 


(2.47) 
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An  expression  for  the  vciometric  emission  coefficient, 
must  be  obtained  in  order  to  evaluate  the  integral 
in  equation  (2.46).  Since  there  are  no  magnetic  fields 
and  hence  no  cyclotron  radiation,  and  since  reccmibination 
of  atoms  is  negligible  at  thermonuclear  temperatures, 
the  only  non-neglible  radiative  process  is  ion  •electron 
bremsstrahlung .  Therefore  the  volumetric  emission  coeffi¬ 
cient  is  just  the  power  due  to  bremsstrahlung  per  unit 
volume  per  unit  s;olid  angle  per  unit  frequency,  which  can 
be  obtained  from  equation  (2.9)  by  changing  variables  from 
wavelength  to  frequency.  If  the  resulting  volumetric  emission 
coefficient  is  integrated  over  all  frequency  and  solid  angle, 
(assuming  isotropic  emission),  one  obtains  an  e;q>ression 

for  the  power  per  unit  volume  due  to  bremsstrahlung  of  all 

28 

frequencies,  (Glasstone  and  Lovberg  ) ; 

=  1.57  X  10"^^  n^  '^(n.  (^K)  ^ -  .  (2.48) 

i  cm  sec 

The  radiation  term  in  equation  (2.43)  is  just  ^  j  P  dx'  . 

In  summary,  the  shock  structure  equations  for  a  radia¬ 
tive  thermonuclear  shock  in  an  optically  thin  atmosphere  are 

/’u  =  A  ,  (2.18) 

p  +  /-  U^  =  B, 


(2.19) 


I 


Pbr«*'  =C, 


(2.49) 


(2.50) 


and  (n^  u)  =  i  =1,  2,...,  N,  (2.51) 

where  A,  B,  and  C  are  defined  by  equations  (2.21),  (2.22), 

and  (2.23)  and  where  the  expression  for  P  is  given  by 

BR 

equation  (2.48).  The  remaining  undefined  quantities,  the 
thermonuclear  energy  release  term,  q^^,  and  the  rate  func¬ 
tions,  W^,  depend  on  the  specific  nuclear  chemistry  and  will 
be  discussed  in  the  next  chapter. 

The  shock  structure  equations  have  now  been  developed 
as  far  as  possible  without  reference  to  a  particular  reac¬ 
ting  gas  mixture.  In  the  next  chapter,  the  nuclear  chemistry 
appropriate  to  three  different  reacting  gas  mixtures,  wilJ 
be  discussed  and  incorporated  into  the  equations. 


i 


C3JAPTER  III 


SHOCK  STRUCTURE  EQUATIONS  FOR  THREE  DIFFERENT 


GASES 


The  general  set  of  radiative  thermonuclear  shock  struc¬ 
ture  equations  derived  in  Chapter  II  will  now  be  developed 
to  apply  to  thermonuclear  shocks  in  (1)  deuterium-tritium 
gas  mixtures  with  equal  numbers  of  deuterium  and  tritium 
ions,  (2)  a  pure  tritium  gas,  and  (3)  hydrogen-tritium  gas 
mixtures.  The  three  cases  will  be  discussed  in  succession. 
After  an  examination  of  the  properties  and  cross  sections 
of  the  nuclear  reactions  occurring  in  the  gas  mixture,  the 
equations  that  govern  reaction  kinetics  will  be  derived. 

The  particular  set  of  shock  structure  equations  for  the  gas 
mixture  will  then  be  put  into  a  convenient  dimensionless 


form. 


3.1  Fusion  Reactions  in  Deuterium-Tritium  Mixtures 

The  following  reactions  occur  in  deuterium-tritium  gas 

39 

mixtures  (Wandel,  Hesselberg  Jensen  and  Kofoed-Hansen  ) 


1.  D  +  D  -^T  +  p, 

2.  D  +  D— »He^  +  n, 

3.  T  +  D— *He^  +  n, 

4.  He^  +  D— »He^  +  p. 
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5  .  T  +  T  +  2n, 

6.  +  He^ — >  He^  +  2p, 

3  4 

7.  T+He  —  'He  +p  +  n. 

Reactions  5  through  7  have  smaller  cross  sections  at  most 
energies  than  reactions  1  through  4  and  are  usually  neglected 
in  the  treatment  of  thermonuclear  problems .  A  graph  of  the 
cross  sections  for  reactions  1  through  4  for  deuteron  energies 
up  bo  1  Mev  is  shown  in  figure  6. 

The  kinetics  of  deuterium  reactions  (reactions  1 
through  4)  are  complicated  because  deuterium  reacts  with 
its  reaction  products,  tritixim  and  helium-3,  Furthermore 
neutrons  and  protons  are  both  intermediate  and  end  products . 
Although  problems  involving  the  full  set  of  deuterium  reac¬ 
tion  kinetics  equations  have  been  solved,  it  is  considerably 
more  difficult  to  deal  with  a  coupled  set  of  fluid  and 
kinetics  ecuations  for  deuterium. 

It  has  been  assumed,  in  the  treatment  of  the  thermo¬ 
nuclear  shock  structure  in  a  deuterium-tritium  mixture  with 
erual  numbers  of  deuterium  and  tritium  ions,  that  only  the 
T(d,n)  He*",  (T-D) ,  reaction  occurs.  This  assumption  is 
well  justified  for  temperatures  below  100  kev  as  reaction 
probabilities  for  T-D  reactions  are  100  times  greater  than 
those  for  the  D-D  reactions  and  10  times  greater  than  that 
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FIG.  6  CROSS  SECTION  FOR  DEUTERIUM  FUSION  REACTIONS 
[glasstone  and  LOVBERG^®] 
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for  the  reaction  He  (d,p)  He  ,  (He  -  D) ,  (see  figure  7). 

For  Mev  temperatures,  the  reaction  probabilities  for  the 
three  reactions  are  comparable.  However,  since  reaction 
rates  depend  both  on  reaction  probcdjilities  and  the  pro¬ 
ducts  of  number  densities  of  the  reacting  nuclei, 

.  »  3 

(W  =  ^1^2  ^  ^  ^  reactions  can  be  neglected  because 

of  the  small  He^  number  densities.  The  ratio  of  the  reaction 
rates  for  T-D  and  D-D  reactions,  assuming  they  have  the  same 
probability  is 


2n„ 

TD  _ T 

DD  D 

The  triton  population  is  enriched  at  the  expense  of  the 

deuteron  population  as  tritons  are  born  in  D-D  reactions, 

so  as  the  reactions  proceed,  the  T-D  reaction  becomes  more 

important  than  the  D-D  reaction. 

The  T-D  reaction  cross  section  has  a  wide  resonance 

whose  peak  is  at  0.107  Mev  laboratory  energy  (see  figure  8) . 

The  lifetime  of  the  associated  state  is  2-3  times  the  natural 

nuclear  lifetime,  the  time  for  a  deuteron  to  traverse  a  tri- 
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tium  nucleus  (Kaplan  ) .  The  main  reaction  mechanism  at 
these  energies  is  formation  of  a  He^  compound  nucleus 


which  decays  into  a  neutron  and  alpha  particle.  At  higher 
energies,  stripping  reactions  dominate;  an  alpha  particle 
and  neutron  are  formed  directly  when  the  deuteron  strips 


I  2  4  6  BK)  OO 

kT(keV) 


1000 


FIG.  7  REACTION  PROBABILITIES  FOR  DEUTERIUM  FUSION 
REACTIONS  BASED  ON  MAXWELLIAN  DISTRIBUTIONS. 
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a  deuteron  from  a  triton  (or  a  triton  strips  a  proton  from 
a  deuteron) . 

The  only  reaction  between  a  deuteron  and  a  triton 

4 

that  is  observed  at  Mev  temperatures  is  T(d,n)  He  .  17.6 

Mev  in  released  in  this  reaction.  The  threshold  energies 

3 

for  two  other  reactions,  T(d, 2n)  He  and  T(d,pn)  T,  are 

2.226  and  2.991  Mev  (center  of  mass)  respectively.  These 

reactions  are  not  observed  for  deuteron  laboratory  energies 

41 

below  about  12-14  Mev  (Ajzenberg-Selove  and  Lauritsen  ) . 

According  to  classical  mechanics,  no  reactions  should 
occur  for  deuteron  energies  less  than  the  Coulomb  barrier  energy 
which  is  ^  660  kev  for  a  tritium  nucleus.  Quantma 

mechanically,  deuterons  with  energies  less  than  the  barrier 
energy  have  finite  probabilities  of  tunneling  through  the 
barrier.  The  probability  of  penetrating  the  barrier,  the 
Gamov  factor,  is 


P  =  exp 


ZjZze^  (^2uf\ 

fi  iEV  / 


(3.1) 


z  z  © 

where  r  ■  1  2  is  the  classical  distance  of  closest 

E  ’  ^ 

approach  in  a  collision  and  \  =  '(2  «<  E)^  Broglie 

wavelength.  The  maximum  reaction  cross  section  for  an 


iu 


s-wave  interaction  (neglecting  Coulomb  repulsion)  is  "n  ^  . 
A  plausible  formula  for  the  reaction  cross  section  for 


V 


/ 
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energies  much  less  than  the  barrier  height,  is  the  product 
of  the  maximum  collision  cross  section  and  the  Gamov  barrier 
penetration  factor.  The  Gamov  reaction  cross  section 


A*  n*  \ 

fT  exp(-  ,  (3.2) 

E  '  (E)" 

has  this  energy  dependence.  Although  the  theoretical  value 
of  B*,  given  in  equation  (3.1)  is  sometimes  used,  the  values 
of  A'  and  B*  are  usually  obtained  from  the  best  available 
low  energy  data.  The  Gamov  cross  section  fits  the  experi¬ 
mental  data  well  for  deuteron  energies  less  t*ian  20  kev. 


4 

3,2  Reaction  Probabilities  for  T(d.n)He  Reactions 

The  reaction  probability  for  T-D  reaction  must  be 

computed  numerically  for  the  temperatures  occurring  in 

thermonuclear  shocks .  Cross  sections  for  energies  between 

1  and  10  kev  were  computed  from  the  Gamov  formula  with 

A'  =  2.19  X  10^  barns  kev  and  B'  =  44.24  (kev)^'^^.  106 

other  values  of  the  cross  section  for  energies  between  10 

kev  and  10  Mev  have  been  taken  from  the  experimental  data 
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shown  in  figure  8  (LA-2014)  .  was  evaluated  with 

a  Simpson's  1/3  rule  adapted  to  unevenly  spaced  data  points, 
(Uneven  spacing  of  data  points  made  it  possible  to  deter¬ 
mine  the  number  of  points  according  to  how  rapidly  the 
cross  section  was  changing  with  energy,)  A  graph  of  <<rv> 


for  kr  between  10  kev  and  3  Mev  is  shown  in  figure  9, 

Reaction  probabilities,  ^  for  several  values  of 

kr  between  1  kev  and  150  kev  were  compared  with  the  results 
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of  earlier  calculations  done  by  Wandel  et  al  ,  Thompson  , 

38 

and  those  quoted  by  Glasstone  and  Lovberg  (p.  18,  based 
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on  a  combination  of  the  work  of  wandel,  Thompson,  and  Tuck  ). 
These  results  are  presented  in  Table  2. 

The  cross  sections  can  be  read  accurately  from  the 
graph  to  two  decimal  places  only.  Therefore  agreement  to 
two  decimal  places  between  our  results  and  those  obtained 
by  the  others  is  all  that  is  meaningful.  The  experimental 
data  used  by  Thompson  and  Glasstone  and  Lovberg  is  not 
the  data  that  was  used  in  our  calculation.  Wandel  et  al, 
however,  used  the  data  presented  in  IA-2014  and  their  results 
agree  with  ours  to  two  decimal  places.  Disagreement  bet¬ 
ween  our  values  of  and  the  values  computed  by  Thompson 

and  Glasstone  and  Lovberg  for  kT=  2  kev  is  explained,  by  the 
fact  that  the  contribution  of  the  cross  section  for  lab 
energies  below  1  kev  was  neglected  in  our  calculation. 
However,  for  the  purposes  of  this  study,  the  value  of  <o‘> 
at  kr  =  2  kev  is  never  used. 


Glasstone  and  /SSL.\ 
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TABLE  2.  Reaction  Probability  for  T(d,n)He  Reaction 


3.3  Kinetics  Equations  for  Deuterium-Tritium  Reactions 

If  only  the  T-D  reaction  occurs,  a  reacting  deuterium- 
triti\iro  mixture  is  made  up  of  tritons,  deuterons,  neutrons, 
alphas  and  electrons.  The  total  pressure  of  the  mixture  is 


p  =  ("t  +  "d  +  "n  n.  +  • 


(3.3) 


The  number  densities  of  the  different  species  can  be  related 
to  each  other  by  the  equations  of  charge  conservation. 


n  =  n_  +  n_  +  2n.  , 

e  T  D  ' 


(3.4) 


electron  „ur*er  flux  conservation  (electrons  do  not  partici- 
P^te  in  the  reactions). 


n  u  =  p. 


(3.5) 


and  mass  conservation. 


^  T  "V  *  "d  "d  "n  ”»  *  )u  -  A.  (3.6) 

P,  the  electron  flow  constant  and  A,  the  mss  flow 

>  wie  mass  flow  constant 

can  be  evaluated  from  the  initial  conditions.  Since  the 

initial  mixture  contains  deuterium  and  tritium  ions  in 
equal  numbers. 


n  =  2n 
e^  . 


(3.7) 


and  therefore 


f  =  2n  u 

T  I  ^ 


(3.8) 


A  =  2. 


3  "t.  "V  ^ 

1  1 


(3.9) 


A  combination  of  eouations 


(3.4)  -  (3.9)  yields 


("t  -  n  )  .  2(„,  .  „ 


(3.10) 


Since  n  =  n  av,.^  - 

Since  It  is 

do  not  diffuse. 


assumed  that  the 


species 


D  T 


(3.11) 


for  all  X.  Therefore,  it  follows  from  equation  (3.10)  that 


N  ^ 


(3.12) 


for  all  X  also. 


The  continuity  equations  for  the  species  are  not  all 


independent.  Tritons  and  deuterons  are  destroyed  at  the 
same  rate,  so  “  ^D*  ^^J^thermore,  neutrons  and  alpha 
particles  are  created  at  the  same  rate,  so  •  Since 

mass  is  conserved  in  the  reactions,  ^  m.  W.  =0,  and  there- 
fore  *=  “^T*  expected,  the  rate  of  triton  and  deuteron 
destruction  is  equal  to  the  rate  of  neutron  and  alpha  creation 
There  is  only  one  independent  continuity  equation  (which  we 


choose  to  be  the  triton  equation) . 

In  order  to  eliminate  number  densities  from  the  flow 
equations,  we  define  a  dimensionless  variable  ,  the  ratio 
of  the  density  of  the  reaction  products  to  the  total  density 
at  the  point  x: 


n„  itL,  +  n  m 
N  N  >4  ' 


(3.13) 


where 


n„  m„  +  n„  iiL  +  n  m 
D  D  N  N  A 


(3.14) 


I 


f 
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Before  the  gas  begins  to  react,  n„  =  n,  =  0,  so  oC=  0. 

After  the  gas  has  reacted  completely,  =  0,  so  (<=  1. 
A  combination  of  equations  (3.11)  -  (3.14)  gives  expressions 
for  the  neutron  and  triton  number  densities: 


n 


N 


3 

5  nvj. 


(3.15) 


and 


n_ 


3  F{-L-<k\  . 

5 


(3.16) 


With  these  expressions  for  n^  and  n^^,  the  state  equation 
(3.3)  becomes, 


P 


12  /:  k_ 
5  n^p 


T  , 


(3.17) 


12k  R 

which  is  independent  of  .  Since  r~  =  T77  >  the  state 

5m,j,  5/4 

equation  can  be  written  in  its  usual  form  for  a  mixture: 


P 


T 


where 


(3.18) 


In  general,  a  state  equation  depends  on  the  degree  of 

reaction  only  if  the  number  of  particles  before  the  reaction 

is  not  equal  to  the  number  of  particles  after  the  reaction. 

3  2  1  4 

For  the  reaction,  T  +  D  — #  n  +  He  ,  the  number  of  particles 

before  and  after  the  reaction  is  the  same;  for  the  reaction, 

3  3  4  1 

T  +  T  — »  He  +  2n  ,  however,  the  particle  number  goes  from 
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two  to  three.  If  the  particle  number  does  not  change  in  the 
reaction,  the  mean  molecular  weight  of  the  reacting  mixture 
remains  constant. 

If  we  rewrite  the  continuity  equation  for  tritons  in 
terms  of  <  ,  using  equation  (3.16),  it  becomes 


2  A  ^ 

5  m  dx 
T 


=  (cm"^  sec“^).  (3,19) 


The  thermonuclear  power  per  unit  volume  term, in  the 
differential  form  of  the  energy  equation,  (equation  (2.17)), 
is  —  (  ^  can  also  be  written  as 

the  product  of  the  reaction  rate,  W,  and  the  energy  released 
per  reaction,  Equating  these  two  expressions  for 

we  see  that 


A 


dx 


dj( 

dx 


(3.20) 


If  both  sides  are  integrated  with  respect  to  x,  subject  to 
the  condition  that  =  0  for  '<=0,  equation  (3.20) 
becomes 


Q  <  , 


Q  ^ 


3 

5 


(3.21) 


where 
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It  is  easily  verified  that  Q  =  — ■ —  ,  the  maximum  energy 

m^  +  lUj 

released  per  unit  mass  of  initial  mixture  if  it  reacts  com¬ 
pletely. 


3 .4  Dimensionless  Form  of  Shock  Structure  Equations  for 
Deuterium-Tritium  Mixtures 

The  shock  structure  equations  can  be  put  into  dimen¬ 
sionless  form,  if  they  are  multiplied  by  an  appropriate 
combination  of  the  mass,  momentum,  and  energy  flow  constants, 
A,  B  and  C,  The  dimensionless  temperature  is  defined  by 


t  3 


k  T  H 


mj 


where 


(3.22) 


(3.23) 


and  m^  is  the  average  ion  mass  for  the  inital  mixture.  For 
the  deuterium-tritium  mixture,  with  equal  numbers  of  deuterium 
and  tritium  ions,  m^  =  (m^j  +  n\j,)/2.  Therefore, 


T  =  T  H 

"(»nj3  +  '\)/2 


R  T  H 
2.5 


(3.24) 


The  dimensionless  speed,  density  and  pressure  are  defined  in 
the  following  way: 
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(3,27) 


In  stinunary,  the  complete  system  of  shock  structure 


equations  in  dimensionless  form  is 


<^.  u;  =  1  ^ 


Tl  +  u.'  =  1, 


5  TT  U)' 


2  -^  +  Y  -  oCHQ  +  A  =  CH, 


TT=  2  , 


4^  _  1  M 

dx  5  nu 


<T\/>  , 


Where 


-J.Iv.  J 

P  =  2.261  X  10“^^  -  f— )  r*. 

BR  2  \ 

VR  / 

T  ^  / 


(3.28) 


(3.29) 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


The  function  A,  can  be  rewritten  in  a  way  that  makes  its 
physical  significance  apparent  by  changing  variables  of 
integration  from  x  to  oC  and  by  making  use  of  equation  (3.32) 
It  then  becomes 


\  BR  ,  t 
=  HO  r —  d< 

J  TN 


(3.35) 


There  are  two  limiting  cases  of  the  shock  structure  equa¬ 
tions,  one  including  reactions  and  no  radiation,  the  other 
including  radiation  but  no  reactions .  The  non-radiarive  limiting 
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case,  (A  -  0  for  all  x),  was  discussed  previously.  The  non¬ 
reactive  limiting  case  is  obtained  if  we  set  ^  “  0 
X  and  o(,(0)=0.  The  radiation  is  assumed  to  occur  downstream 
from  the  gas  dynamic  shock,  so  the  physical  validity  of  this 
case  is  questionable.  However  the  pure  radiative  solution 
of  the  equations  is  useful  in  providing  a  bounding  curve 
for  the  solution  of  the  physically  more  valid  case  that 
includes  both  reactions  and  radiation. 


3.5  Dimensionless  Parameters  in  the  Energy  Ecmation 

There  are  three  dimensionless  parameters  in  the  energy 
equation,  CH,  HQ,  and  A. 

The  parameter,  CH,  combines  the  constants  of  integration 
of  the  mass,  momentum  and  energy  equations,  it  depends  only 
on  the  initial  Hach  number  of  the  gas,  and  not  on  the 


reaction  chemistry,  i.e.. 


/5  n  “I 

CH  =  -  -  +  - 

U  />  •> 


=  i  .  -5-  ] 

2  V  / 


/Li 


,  2 


1+1 


(3.36) 


The  initial  Mach  numbers,  for  thermonuclear  shocks  are 
>  >  1,  so  CH  =  1/2.  The  value  of  CH  is  therefore  independent 
of  the  initial  conditions. 
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The  parameter,  HQ,  combines  the  mass  and  momentum  flow 
constants  with  the  energy  per  unit  mass  released  in  the 
reactions,  i.e.. 


y  2 

For  very  large  initial  Mach  numbers,  HQ  reduces  to  Q/u^^; 
since  u^^  and  Q  are  related  by  the  Chapman-Jouguet  condition, 
then 

].  9 

HO  =  - - —  =  —  .  (3.38) 

2(  i  -1)  32 

HQ  is  therefore  independent  of  the  reaction  under  considera¬ 
tion  for  shocks  propagating  with  the  C-J  speed. 

The  third  dimensionless  parameter,  A,  is  a  measure  of 
the  relative  importance  of  power  lost  in  bremsstrahlung 
and  power  gained  in  thermonuclear  reactions.  When  A  is 
much  smaller  than  other  terms  in  the  energy  equation,  (A«<1), 
radiative  losses  are  insignificant.  When  the  value  of  A 
approaches  1,  radiative  power  losses  become  comparable  to 
the  power  gained  in  reactions.  Unlike  the  parameters  CH  and 
HQ,  A  is  neither  constant  nor  independent  of  the  reaction 
chemistry. 

3.6  Fusion  Reactions  in  Tritium 


Two  fusion  reactions  can  occur  in  a  pure  tritium  gas. 


1 


4 

T  +  T  — >  2n  +  He  j 
and  2.  T  +  T  — >He^+  n;  He^-- >  He^  +  n. 

11.33  Mev  is  released  in  each  reaction.  If  an  incident 
triton  strips  a  proton  off  a  target  triton,  reaction  1 

4 

occurs,  and  He  and  two  neutrons  are  formed  directly.  The 
neutron  energy  spectrum  is  continuous  since  there  are  three 
bodies  after  the  reaction.  If  an  incident  triton  strips  a 
deuteron  off  a  target  triton,  the  reaction  occurs  in  two 
stages  (reaction  2  ) .  The  neutron  energy  spectrum  is  dis¬ 
crete  because  only  two  bodies  come  off  in  each  stage  of 
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the  reaction  (A jzenberg-Selove  and  I^auritsen  ) . 

T-T  cross  sections  have  been  measured  for  incident  lab 

40 

energies  between  60  kev  and  1.14  Mev  (Govorov  et  al  ).  The 

R 

cross  section  as  a  function  of  energy  may  be  represented  by 
the  function 

7  =  (a  +  b<nE  (kev))  x  10  cn^  (3.39) 

where  a  =  -91.2  +  2.5, 

b  =  (55.8  +  D/inlO  • 

r  increases  monotonically  from  10  mb  at  60  kev  to  82  mb 
at  1.14  Mev  (see  figure  11). 

Below  these  energies  we  can  assume  a  cross  section  of 
the  Gamov  type.  There  is  seme  question  about  the  applica¬ 
bility  of  Gamov  cross  section  to  stripping  reactions.  How¬ 
ever  it  is  the  only  available  cross  sections  for  these  energies 


Y-<- 


and  has  been  used  by  Thompson  .  The  theoretical  value  o£ 

1/2 

the  constant  B'  is  54  (kev)  ,  The  value  of  A'  which  gives 
Govorov's  measured  value  of  cr  at  60  kev  is  * 


A'  =0.6  exp/  —  1 
1(60)7 


barns  kev 


(3.40) 


The  T-T  cross  sections  for  E  4  60  kev  are  shown  in  figure  10 

Jj 


3.7  Reaction  Probabilities  for  T(  t.2n)  He  Reactions 


An  expression  for  the  reaction  probability  for  T-T 

reactions  can  be  obtained  analytically  for  the  temperatures 

that  occur  in  thermonuclear  shocks  (kT  ^  1  Mev) .  In  order 

to  compute  an  assumption  must  be  made  about  the  behavior 

of  the  cross  section  at  energies  for  which  no  experimental 

values  are  available  (E_  <,  60  kev  and  E_  >  1.14  Mev) .  For 

ij  L 

kT  1  Mev,  collisions  between  particles  with  lab  energies 
less  than  60  kev  do  not  contribute  significantly  to 
Therefore,  any  convenient  assumption  about  the  behavior  of 
the  cross  section  at:  these  energies,  can  be  made.  The 
Govorov  cross  section  has  been  extrapolated  to  its  crossing 
point  on  the  energy  axis,  E^  =  exp(-  ^  )•  It  is  assumed 
to  be  zero  for  energies  below  E^.  Collisions  between  parti¬ 
cles  with  energies  above  1.14  Mev  make  important  contributions 
to  <5''>  .  The  Govorov  cross  section  has  been  extrapolated 


*  An  alternative  way  to  match  the  Gamov  and  Govorov  cross 
sections  at  60  kev  would  be  to  choose  values  of  A*  and  B* 
such  that  the  cross  section  and  its  first  derivative  were 
continuous  at  60  kev. 


The  integral  can  be  evaluated  most  simply  by  changing 

variables,  from  E  to  the  dimensionless  variable,  ^  ,  where 

cm 


(3.42) 


Then 

*27  *  r"  ^ 

{Tv’>  « r(a'+bfnfe)e  td€  ,  (3.43) 

where  a'«  a  +  b  Jin(2kT) 

When  the  integral  is  performed,  the  resulting  expression 


for  <  1  v  >  is 


(3.44) 


Written  in  terms  of  our  dimensionless  variables,  becomes 

(3.45) 


L  .  i 

- i  b 

r  2'^  )j 

Since  H  - - ^ -  , 

2(r-l)Q 


E  H 
O 


is  a  dimensionless 


pareuneter  that  depends  on  the  Q  of  the  reaction  and  the  cut¬ 


off  energy  for  the  Govorov  cross  section,  E^. 


kT  ~  100  key 

If  kT  ~  100  kev,  the  contribution  of  collisions  of 
particles  with  lab  energies  below  60  kev  becomes  important. 
The  Gamov  cross  section  must  be  used  for  energies  below  60 
kev.  There  are  now  two  contributions  to  The  Gamov 

cross  section  is  integrated  from  some  energy  (which  is  much 
less  than  kT)  to  E  =  60  kev.  The  Govorov  cross  section  is 

Ij 

integrated  from  E_  =  60  kev  to  E  »  a?  ,  The  Geunov  cross 

L  L 

section  was  evaluated  at  37  values  of  lab  energy  between 
0.1  and  60  kev,  and  integrated  over  these  energies  with  a 


Simpson's  rule  for  unevenly  spaced  pivotal  points.  The 
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Integral  over  the  Govorov  cross  section  was  performed 
analytically  with  the  result 

=  .  -  OcT)”'  r<r  (1+  )  exp  ( — ) 

L  M  kT 

where  <t-  =  a  +  b  Jl  n60  and  E  =30  kev. 

o  oo 

The  results  of  the  calculations  for  <Tv'>^  and 

'Gam  '  Gov 

are  shown  in  Table  3,  For  kT  =  1  kev,  the  Gamov  contri¬ 
bution  dominates;  for  kT  =  10  kev,  the  Govorov  contribution 
is  roughly  twice  the  Gamov  contribution;  for  100  kev,  the 
Govorov  contribution  dominates. 


kT (kev) 

Gov 

3 

.  gn  . 

OV>  t  sec^ 

1 

3.56  X  lO"^’^ 

1.04  X  10”^° 

3.56  X  lO"^^ 

5 

-20 

6.64  X  10 

1.88  X  10”^° 

-20 

8.52  X  10 

10 

-19 

2.11  X  10 

-19 

3.99  X  10 

-19 

6.10  X  10 

20 

-3*' 

2.44  X  10 

2.23  X  10“^® 

2.47  X  10“^® 

50 

1.29  X  lO"^^ 

8.97  X  10“^® 

9.10  X  10“^® 

100 

5.84  X  10”^° 

-17 

1.93  X  10 

-17 

1.93  X  10 

200 

-20 

2.34  X  10 

3.68  X  lO"^^ 

3.68  X  lO"^^ 

500 

6.38  X  lO”^^ 

7.83  X  10“^^ 

7.83  X  lO"^^ 

900 

2.70  X  lO"^^ 

_ 

1.22  X  10“^^ 

1.22  X  lO"^^ 

TABLE  3.  Contributions  of  Gamov  and  Govorov  Reaction 
Probabilities  to 


An  approximate  formula  for  <<\J>  for  kT  -4  10  kev,  is 


44 

(Thompson  ) 


n  o  -2/3 

<<rM>  =  0^8-X.  10,,  AB„  .  „ 

app 


-3/21^ 

\^4W/ 


2\  1/3 


cari/sec,  (3.47) 


where,  for  T-T  reactions. 


=  the  reduced  mass  in  a.m.u. 
A 


Thompson  has  used  the  Gamov  constants.  A*  =59.9  barns  kev 

*1/2 

and  B'  =  668  (kev)  .  However,  the  values  of  for 

kT  =  1  kev  and  kT  =  10  kev,  given  in  his  paper,  are  incon¬ 
sistent  with  the  values  of  obtained  with  his 

app 

values  of  the  Gamov  constants.  The  Gamov  constants, 

1/2 

A'  =  54  bams  kev  and  B'  =  630  (kev)  ,  have  been  used  in 
this  study.  The  values  of  <‘^V>  are  compared  with  the 
results  of  the  numerical  integration  of  the  Gamov  cross 
section  for  kT  =  1  kev  and  10  kev  in  Table  4.  The  agree¬ 
ment  between  the  numerical  results  and  the  analytical 
results  is  quite  good  since  the  analytical  formula  is  only 
approximate.  This  constitutes  a  check  on  the  numerical 
integration  for  small  values  of  kT. 


y  A  A'=  54 

^^^Gam  B ' »  630 

A'=  54 

<'^>app  B'=  630 

3.56  X  10"^^ 

6.10  X  10"^^ 

-23 

3.32  X  10 

7.00  X  10“^^ 

A'=59. 
^v)>app  B'=668 


0.795  X  10“^^ 
3.69  X  10“^^ 


TABLE  4.  Reaction  Probabilities  for  Low  Temperatures 

The  T-T  reaction  probability  for  kT  between  12  kev  and 
2  Mev,  is  shown  in  figure  9,  along  with  the  T-D  reaction 
probability.  The  T-D  reaction  probability  rises  to  a  maxi¬ 
mum  value  at  70  kev  and  then  decreases  with  increasing 
temperature.  This  shape  shows  the  influence  of  the  resonance 
in  the  T-D  reaction  cross  section.  The  T-T  reaction  probabil¬ 
ity,  however,  rises  inonotonically  with  increasing  temperature 
Although  T-D  reaction  probabilities  are  about  100  times  larger 
than  T-T  reaction  probabilities  for  kT  10  kev,  the  two 
become  comparable  at  —  2  Mev. 


3.8  Kinetics 


or  Tritium  Reactions 


A  reacting  tritium  gas  mixture  consists  of  electrons, 

5 

tritons,  alfdias  and  neutrons.  He  has  a  lifetime  of  the 
-21 

order  of  10  sec  so  it  need  not  be  included  as  one  of  the 
species  of  the  mixture.  The  total  gas  pressure  is 


P  =  +  Ik,  +  n  +  n^)kT. 


(3.49) 
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The  ntuiiber  densities  of  the  species  are  not  independent. 
Since  the  gas  has  local  neutrality. 


n^  =  n^  +  2n^  .  (3.50) 

No  electrons  are  created  or  lost  in  the  reactions  so  the 
electron  number  flux  is  conserved: 
n  u  =  F, 

where  P  =  n  u,  .  (3.51) 

Ti  1 


The  total  mass  flux  is  also  a  constant, 

+  "n  ”*N  n^  )  u  =  A, 


(3.52) 


where 

Since  A 


A  =  n_ 


l^^l  * 


n^F,  it  follows  that 


(3.53) 


The  electron  number  density  profile  has  the  seune  shape  as 
the  total  density  profile  (which  is  reasonable  since  the 
electrons  are  inert) .  The  relationship  between  the  alpha 


and  neutron  number  densities. 


n. 


n. 


N 


(3.54) 


follows  when  the  equations  of  electron  flux  conservation 
and  mass  flux  conservation  are  combined  with  the  condition 
for  local  charge  neutrality. 

The  continuity  equation  for  tritons  is 


-  ,^..,.jKjF'»;>J-,4w<»^-fS  ■■'  E ' 
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|-  (”t  .  Mj  =.J^  ct^>  . 


(3.55) 


W^l  is  twice  the  reaction  rate  for  tritons  since  two 
tritons  are  lost  in  each  reaction.  The  continuity  equation 
for  alpha  particles  is 


^  <"*  “>  =  . 
dx 


(3.56) 


The  rate  of  neutron  creation  is  twice  the  rate  of  alpha 
particle  creation  so 

Wjj  =  2W^  .  (3.57) 


Since  the  total  mass  flux  is  conserved.  Zm.W.  =0.  and 

ill 

therefore 


T  N 


(3.58) 


It  is  convenient  to  eliminate  number  densities  from 
the  continuity  and  state  equations.  c<.  is  defined  in  the 
usual  way  as  the  ratio  of  the  density  of  the  reactions  pro¬ 
ducts  to  the  total  density. 


<  = 


_  n^  m^ 


Z’ 


(3.59) 


where  f  -  \  *^1  *  (3.60) 

If  equation  (3.59)  is  solved  for  n^^  and  n^,  the  resulting 


expressions  are, 

"n 


(3.61) 


The  state  equation  can  be  rewritten  in  terms  of  o<^  .  With 
Ic  R 

the  relation  —  "  T  ^  it  becomes. 


p  =fR  T  (2+  1). 
3  2 


(3.63) 


This  equation  is  recognized  as  the  familiar  equation  for  a 


mixture  of 


gases,  p  =  P  ^  T,  if  it  is  noted  that  the 


mean 


molecular  weight  per  particle  of  the  mixture,  W,  instead 


of  being  constant,  varies  with  oC  ,  i.e,, 
W  =  3/(2  +  o</2) . 


(3.64) 


The  set  of  dimensionless  variables  that  were  used  to 
make  the  deuterivun-tritivun  structure  equations  dimensionless, 
are  now  used  to  make  the  tritium  structure  equations  dimen¬ 
sionless.  The  appropriate  dimensionless  temperature  is. 


T  = 


k  T  H  _  R  T  H 
m_  ~  3 


(3.65) 


The  mass  and  rooroentvun  equations  are  the  s^uae  for  the  deuterivun- 
tritium  and  pure  tritium  gases.  The  specific  reaction  chemistry 


is  contained  in  the  state  equation, 

■n  =  (2+  ~)  , 

the  continuity  equation, 

dx  itL,  U'2  ^  '  » 


(3.66) 


(3.67) 


(  .j. ..  -»»»»« ■*- 
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and  thtt  radiation  term  in  the  energy  equation, 

< 


d  -  HQ  J  fgjR  d  oc' 


Where 


and 


0  ^TS 


P__  «  1.57  X  10 
oR 


P_-  «  AQ  d< 

™  at 


> 

(3.68) 

A  B  H 

2 

(if  r* 

W  w* 

,  (3.69) 

(3.70) 

4 

3,9  Kinetics  Equations  for  T(t.2n)He  Reactions  in  Hvdroaen- 
Tritium  Mixtures 

1£  hydrogen  is  added  to  a  tritium  gas,  the  Chapman- 
Jouguet  speed  is  reduced.  An  exaimination  of  the  shocX 
structure  equations  for  a  mixture  of  hydrogen  and  tritium 
if  of  interest  in  connection  with  shock  tube  experiments. 

The  C-J  speed  for  such  mixtures  can  be  reduced  to  a  value 
which  may  be  obtainable  in  an  electromagnetically  driven 
shock  tube. 

Let  6  be  a  parameter  that  measures  the  initial  density 
ratio  of  protons  to  tritons: 

^  i  •  (3.71) 

Since  reactions  probabilities  for  proton  fusion  reactions  are 

small ,  the  only  energy  release  comes  from  the  tritium  reactions 

Q  is  the  energy  released  per  reaction,  the  energy  released 
R 

per  gram  of  gas  of  the  initial  mixture  is 

I 


'■fy,/. 


Q  »  ?R 


(3.72) 


1+  e 


The  gas  consists  of  five  species t  electrons,  protons, 
tritons,  alphas,  and  neutrons.  The  species  nuinber  densities 
can  be  related  by  charge  conservation. 


=  ”h  -v  + 


conservation  of  electron  flux. 


n^  u  -  F, 


conservation  of  proton  flux. 


"h"h  “  ■  °  ' 

and  conservation  of  mass. 


(3.73) 


(3.74) 


(3.75) 


where  F,  A,  and  G  are  constants  of  integration  which  can  be 
evaluated  from  the  initial  conditions.  The  initial  electron 


density  is 


\  \  '\  >• 


(3.76) 


Therefore, 


P  =  n^^  (1+3  6 )  u^. 


®  \  “l  \“l  * 


(3.77) 


(3.78) 


and  ^  "Si  '^l  "  ®T  ^  ^'^1* 

Only  one  of  the  three  constants  of  integration  is  indepen¬ 


dent,  i.e.. 


(3.80) 


•  H  ‘ 
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and 


i  A. 

(1+t)  «\j. 


(3.81) 


If  w  solve  for  1/u  in  the  above  equations,  %re  obtain 

1  .£  .2a. 

u  A  P  G  * 


(3.82) 


This  relationship  can  be  used  to  get  two  different  esqpressiona 
for  ng. 

^  -  n,  +  3  t  (n_  +  + 


and 


”h  2  ~  '"T  2 

lig  -  6  (3n^  +  lijj  +  4n^  ) . 


*V  )> 


(3.83) 

(3.84) 


Setting  the  two  expressions  for  n^  equal  to  each  other, 
we  see  that 


2a. 


^(n.  -^) 


(3.85) 


Since  6  is  an  arbitrary  input  parameter,  equation  (3.85) 
holds,  (as  in  the  pure  tritium  case,  to  which  it  reduces 
vdien  t  =0) ,  only  if 


2 


n 


U.  » 


(3.86) 


for  all  X. 


The  total  pressure  of  the  gas  mixture  is 

p»  |^2(l+3  fe  )  (n^  +  rijj)  +  ^  J  kT  .  (3.87) 

It  is  convenient  to  eliminate  number  densities  from  the 
state  and  continuity  equations  and  express  them  in  terms  of 
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i 


'  4 


r 


V.. 


some  appropriate  reaction  parameter  o(  .  One  such  choice  of 
o(  puts  the  state  equation  into  the  same  form  as  for  the 
pure  tritium  case.  Let 


<  =  (“»•»;  . , 


(3.88) 


.  ^  -  1+ 


An  effective  triton  mass. 


(3.89) 


(3.90) 


(3.91) 


can  be  defined.  Expressions  for  the  neutron  and  triton 


number  densities  follow  from  equations (3 .88)  and  (3.90): 


(3.92) 


1  - 
1+3  fe 


(3.93) 


The  sum  of  the  neutron  and  triton  number  densities  is  inde¬ 
pendent  of  o(.  .  Ihe  number  densities  of  the  two  inert  compo¬ 
nents  (protons  and  electrons)  of  the  gas  mixture  are  propor¬ 


tional  to  (n^  +  Ujj),  i.e., 

Ujj  »  3  t(n^  +  njj) , 

and  -  (1+3  )  (.n^  +  n^) 


(3.94) 


(3.95) 


Xf  equations  (3.92)  end  (3.93)  are  substituted  into  the 


state  equation,  it  becomes 

p  [2**  ]  . 


(3.96) 


This  is  identical  to  the  state  equation  for  pure  tritium 
if  we  replace  the  triton  mass  by  its  effective  mass.  It 
must  be  noted  that  vdien  the  reactions  are  complete, 
n^  =  0  and  therefore  from  equation  (3.90),  d.  *  • 

The  energy  equation  for  the  hydrogen-tritium  mixture 
is  the  same  as  equation  (3.30).  The  expression  for  A  must 
be  modified  to  contain  the  appropriate  reaction  chemistry, 
and  the  constant  Q  in  the  thermonuclear  energy  release  term 
must  be  redefined.  The  thermonuclear  energy  release  term, 
oC  Q,  has  the  value  Q  for  a  tritium  gas  when  all  the  reac¬ 
tions  are  completed.  The  maximum  value  of  oC  for  the  hydro- 
gen-tritium  mixture  is  maximum  energy  released 

pe^  unit  mass  of  the  initial  mixture,  (given  by  equation 

Or  1 

(3.72)),  is  Q  ■- —  ; " ■ ;  .  We  define  a  constant  Q 

*  2iaji  1+  €  erf 

such  that  Q  --  approaches  Q  as  ci  approaches  ..JL.  . 

It  follows  that 


°eff  “ 


(3.97) 


The  flow  variables  are  made  dimensionless  in  the  same 
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way  as  those  for  the  pure  tritium  case.  The  dimensionless 
temperature  becomes 


-  h  T  H 


(3.98) 


As  in  the  state  equation^  the  mass  of  the  triton  has  been 
replaced  by  its  effective  mass. 

The  dimensionless  mass^  state  and  momentum  equations 
for  the  hydrogen-tritivim  mixture  are  the  same  as  those  for 
pure  tritium.  The  dimensionless  energy  equation  is 
2 

f  f  =  CH, 

where  L  is  defined  in  the  usual  way  and 


P 

BR 


-27  /  3t  7 

1.57  X  10  '  ^  I  Rf 

^  V  R  / 


i: "" 


(3.100) 


The  continuity  equation  is. 


M  M 
dx  in^ 


<^\J> 

X 

U) 
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(3.101) 
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CHAPTER  IV 

SHOCK  STRUCTURE  CURVES  IN  PHASE  SPACE 

4.1  Non-Radiative  Shock  Structure  Curves  in  Phase  Space 

The  system  of  non-radiative  shock  structure  equations, 
developed  In  Chapter  III  for  three  different  gas  mixtures, 
includes  four  algebraic  equations,  (mass,  momentum,  energy 
and  state  equations)  and  one  non-linear,  first  order,  dif¬ 
ferential  equation,  (the  continuity  equation) .  In  order 
to  integrate  the  differential  equation  for  the  shock  pro¬ 
files  in  real  space,  all  the  flow  variables  except  one 
must  be  eliminated.  The  algebraic  equations  are  solved 
for  relationships  between  the  variables.  These  relation¬ 
ships  are  called  the  phase  space  shock  curves. 

A  relationship  between  the  dimensionless  speed,  OJ  , 
and  the  degree  of  fusion,  oC  ,  can  be  obtained  from  the  non- 
radiative  mass,  momentum,  and  energy  equations.  A  combina¬ 
tion  of  equations  (3.28)  and  (3.29)  gives  an  expression  for 
—  ,  which  can  be  substituted  into  equation  (3.30).  The 
resulting  equation  is 

Vx.)  _  lu)  +  "I  (CH  +  HQ  )  =  0,  (4.1) 


or 


<  --i-  /^-ch|  =  • 

^  HQ  i  32  y'  HQ  V  8  / 


(4.2) 


Since  CH  =  1/2  and  HQ  *=  9/32,  equation  (4.2)  becomes, 

'<•-  1  =  -  ^  (w-  ij^  .  (4.3) 

Equation  (4.3)  is  independent  of  the  specific  reaction 
chemistry  since  the  dimensionless  mass,  momentum  and  energy 
equations  are  identical  for  all  three  cases  considered  in 
Chapter  III .  The  u.’  -  cC  curve  is  a  parabola  that  opens  to 
the  left  with  an  axis  of  symmetry  parallel  to  the  oC-axis. 

The  vertex  of  the  parabola  is  located  at  (uj  ,  oC  )  =(l ,  1). 

O 

The  u)  -  oC  relation  is  plotted  in  figure  12.  (Points  1, 

2  and  3  are  the  same  as  in  figure  2.) 

The  two  branches  of  the  lo  -  curve  join  at  the 
vertex  of  the  parabola.  The  lower  branch,  (  )  ,*  corresponds 

physically  to  a  gas  dynamic  shock  followed  by  a  C-j  deflagra¬ 
tion.  It  is  the  fcuniliar  solution  assumed  by  the  von  Neumann 
model.  Ahead  of  the  shock,  lO  is  equal  to  1.  Behind  the  ■ 
shock,  U)  is  equal  to  1/4,  while  behind  the  deflagration, 
has  the  value  5/8.  The  upper  branch,  (u.\),  is  also  a 
solution  of  the  non-viscous  non-thermal  conducting  shock 
structure  equations.  It  corresponds  physically  to  a  C-J 
detonation.  is  a  continuous  function  of  ,  monoto- 

nically  decreasing  as  nC  increases  from  0  to  1.  The  C-J 
detonation  has  not  been  observed  experimentally  and  is 
therefore  not  considered  in  most  theoretical  treatments. 
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FIG.  12  DIMENSIONLESS  SPEED,  w,  VS.  DEGREE  OF  FUSION 
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26 

(e .g . ,  von  Karman  ) .  It  is  presented  in  the  j^ase  plane 
diagrams  for  completeness  sake;  shock  structure  curves  in 
real  space  have  not  been  computed  for  this  branch. 

A  relationship  between  the  dimensionless  temperature, 

'X.  >  and  the  degree  of  fusion,  o<  ,  is  obtained  by  com¬ 
bining  the  mass,  momentum,  and  energy  equations  with  the 
state  equation.  The  mass,  momentxun  and  energy  equations 
are  first  combined  to  give  a  relation  between  and  ^ 
for  the  deuterium-tritium  mixture: 

^0^  -  u;  +  2^  =0;  (4.4) 

and  a  relationship  between  KjJ  ,  X  and  oC  for  the  pure 
tritium  gas: 

-  U.'  +  ^2+  ijir  =  0.  (4.5) 

Combining  these  relations  with  the  appropriate  energy  equa¬ 
tion,  we  obtain  a  relationship  between  X.  and  oC  .  The  X  -oC 
relation,  unlike  the  W-<3<  relation,  depends  on  the  reac¬ 
tion  chemistry  because  it  contains  the  state  equation. 

The  relationship  between  temperature  and  degree  of 
fusion  for  the  deuterium-tritium  mixture  is 

2 

16T^  +  +  2  D'  -  aHQAjT+  +  HQot(HQ«C-  |')  =0,  (4.6) 

where  D'  =  1  -  4CH  *»  -  1,  (4.7) 


x  as  a  function  of  «<  is  plotted  in  figure  13.  The 
jump  conditions  for  temperature  are  obtained  by  setting 


I 


«0  and  1  in  equation  (4.6).  When  c(.>0, 

•  2 

+(  I  +  2D»)T  +  "  0* 

and  =  TT,  0. 


(4.8) 


(4.9) 


The  zero  root  comes  from  assuming  that  D*  >  -1.  Actually 
D*  is  not  exactly  -1  and  When  its  correct  value  is  used,  the 
zero  root  of  equation  (4.8)  becomes  ^  10**^.  The  other 
root  is  only  different  from  3/32  in  2d>out  the  sixth  decimal 
place.  When  =  1, 

2 

16r^+(2D'+  i  -  8HQ)r  +  +  HQ  (HQ-  ~  )  =0,  (4.10) 


la] 


and  r+  =  ^  0.117  . 


(4.11) 


The  branch,  ,  (like  corresponds  to  a  C-J  detonation. 

The  other  branch,  corresponds  to  a  gas  dyneunic  shock 
followed  by  a  C-J  deflagration.  The  two  branches  meet  at  the 
downstream  point. 

■t' -  OK  curve  is  not  any  simple  shape  like  the  «<• 

curve.  However,  the  existence,  location,  and  nature  of  extre¬ 
mal  points  on  it  can  be  determined  analytically,  if  equation 

(4.4)  is  differentiated  with  respect  to  ,  a  relationship 
dr  ,  du> 

between  ^  and  ~  is  obtained: 


_  df  %  du; 

2  —  .  (1.2U.  )  —  , 


(4.12) 


If  ~  =  0,  then  either  ~  =0  or  '*'=  There  are  no  points 
on  the  ^  curve  where  ~  =  0,  but  cu  does  go  through 


I  t\ 


th«  value  1/2  on  the  curve.  Therefore,  there  is  some 
point  on  the  curve  where  ~  «=0 .  If  equation  (4 . 4)  is 

differentiated  a  second  tine,  the  extremum  is  found  to  be 

du)  1 

a  maximum;  (since  T"  0  for  =  t  and  therefore  T^x<0) . 

It  can  be  shown  from  equations  (4.3)  and  (4.4)  that  the 

1  8 

coordinates  of  the  maximum  are  (?',*‘)  =  (’r,  t)  = 

u  9 

(0.125,  0.889). 

The  value  of  the  Mach  number  at  any  point  on  the  « 

curve  is  given  by  the  easily  verifiable  reli'ition 

vP"  2 

=  2  JT  M  .  (4.13) 

2  13 

The  temperature  maximum  occurs  when  M  =  y  “  J  •  Therefore 

1  2  3 

the  temperature  increases  with  increasing  for  <  ■r 

5  5 

3  2 

and  decreases  with  increasing  for  M  ^  When 

2  3 

M  =  —  ,  the  static  pressure,  p,  is  equal  to  the  kinetic 
2 

pressure,  ^u  .  The  temperature  therefore  rises  with  decrees- 

2 

ing  pressure  when  p  > u  ,  and  falls  with  decreasing  pressure 
2 

when  p4  u  .  The  competing  physical  processes  are  the 
reactions,  which  add  energy  to  the  gas,  thereby  raising  the 
temperature,  and  the  transformation  of  random  thermal  energy 
into  directed  kinetic  energy,  which  lowers  the  temperature. 

As  the  reactions  go  to  completion,  the  reaction  rates  decrease 
and  eventually  the  second  effect  dominates  the  first. 
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The  relation  between  the  dimensionless  temperature,  't  , 
and  the  degree  of  fusion,  <*s  ,  for  pure  tritium  and  hydrogen- 
tritium  mixtures  is 

2 

(4+ ^  r  (4+ < )  2- (  ^  -2HQoC  )  + 


+  HQc<  (HQ«t  -  “  )  =  0. 


(4.14) 


7“  as  a  function  of  o<  is  depicted  in  figure  14.  The  Toi 
cross  terms  in  equation  (4.14)  are  of  higher  degree  than  those 
for  the  deuterivuA' tritium  case  because  the  state  equation 
depends  on  oc.  .  The  upstream  (o<=  0)  roots  of  this  equation, 
which  are  the  pre-shock  and  post-shock  dimensionless  tempera¬ 
ture,  are  identical  to  those  for  the  deuterium-tritium  re3.ation 
between  T  and  ot  ,  (equation  (4.6)) .  This  is  reasonable 
since  the  structure  equations  are  made  dimensionless  in  a 
way  that  eliminates  their  dependence  on  the  mean  molecular 
weight  of  the  initial  gas  mixture,  and  since  no  reactions 
take  place  in  the  gas  dynamic  shock.  The  downstream  («<  »  1) 
temperature  satisfies  the  equation 

2 

25r^  +  5(  +  |-  -  2HQ)  r  +  +HQ(HQ-  ^  (4.15) 

vdiich  has  the  double  root. 


m  I 


(4.16) 


For  the  discontinuous  solution  , '?;(  «  =  0)  =7^(oL  =1).  In 
other  words,  the  decrease  in  mean  molecular  weight  per 
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FIG.  14  DIMENSIONLESS  TEMPERATURE  ,  T  ,  vs.  DEGREE  OF 
FUSION  ,  a  ,  FOR  A  DETONATION  IN  TRITIUM 


particle  just  cancels  any  net  rise  in  temperature  across  the 
deflagration  region. 

Hie  existence  of  extremal  points  on  the  T-u  curve 
can  be  shown  also.  If  equation  (4  .5)  is  differentiated 
with  respect  to  o(  ,  then. 


(1-2-)  1^-  (2  -f)  f . 


(4.17) 


This  e^qiression,  unlike  its  counterpart  for  the  deuterium- 

tritium  mixture,  does  not  give  a  simple  condition  for  the 

existence  of  an  extremum  on  the  T  o(  cuirve  •  However,  if 

the  T  -  curve  is  differentiated  directly,  and  we  set 
dr 

~  =0  in  the  resulting  equation,  we  find  an  expression  which 
t'  and  satisfy  at  the  extremal  point.  If  this  expression 
is  substituted  into  equation  (4.14),  the  coordinates  of  the 
extremum  are  found  to  be  »  (0.1036,  0.72),  whicdi  corres¬ 

ponds  to  a  dimensionless  speed,  ou  m  0.427.  It  can  be  shown, 
although  the  algebra  is  more  complex  than  for  the  deuterium- 
tritium  example,  that  this  extremal  point  is  a  maximum.  There 
is  a  second  extremal  point  which  has  been  discarded  because  it 


does  not  correspond  to  a  ^ysically  allowed  value  of  o<>  . 


ve  Shock  Structure  Curves  in  Phase  Space 


When  radiation  energy  loss  is  considered,  the  energy 
equation  changes  frcmi  an  algebraic  equation  to  an  integral 
equation  (see  equations  (3.30)  and  (3.33)).  The  radiative 
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•nergy  loss  term,  A,  can  be  written  as  a  function  of  c< 

(see  equation  3.35).  The  mass,  momentum,  and  state  equa¬ 
tions  can  be  used  to  eliminate  all  variables  except  ^  and 
from  the  energy  equation.  The  resulting  equation, 
relating  ^  and  oC  for  the  radiative  case,  is  the  same 
as  that  for  the  non-radiative  case,  (equation  (4.6)  for 
the  deuterium-tritium  mixture  and  equation  (4.14)  for  the 
tritium  gas),  if  D'  is  replaced  by 


where 


and 


D  =  D'  +  4A  (rfv )  , 


A  =  HQ  R  do( 
J  p 

O 

P  Pth 


(4.18) 


(4.19) 


It  follows  from  equation  (3.32)  and  (3.34)  that  R  can  be 

P 

written  as  the  product  of  a  function  of  (<.  ,  and  a  fum'tion 
of  'X.  ,  i.e.. 


Rp  =  K  A  (o(  )  8(1:  ),  K>0,  (4.20) 

where  A(c<  )  =(1^”  , 

X 

and  B(‘i:  )=  T 

<<r\,')(t) 

An  iterative  method  is  used  to  solve  the  integral  equa¬ 
tion  relating  the  dimensionless  temperature  to  the  degree  of 
fusion.  The  non-radiative  solution  (A  =  0,  for  all  )  is 
assumed  for  the  first  iteration.  The  first  non-zero 
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approximation  to  the  radiative  loss  term,  {c<  ) ,  can  be 
evaluated  with  the  non-radiative  'C'-oC  relation  (if  <<!'/>  is 
a  known  function  of  “T  ) .  A  first  approximation  to  the 
radiative  relation  between  ^  and  o<.  is  then  obtained.  It 
is  used  to  evaluate  A^  (°<.)*  This  iterative  process  is 
continued  until  ^'(o()  does  not  change  to  some  required 
degree  of  accuracy  from  one  iteration  to  the  next.  The 
numerical  procedure  converges  to  three  significant  figures 
in  about  six  iterations.  The  results  of  the  integration 
for  T"  (oC )  are  shown  in  figure  13  for  the  deuterium-tritium 
case,  and  in  figure  14  for  the  pure  tritium  case.  A  (oC ) 
for  both  cases  is  plotted  in  figure  15. 

The  radiative  relation  between  the  dimensionless  speed 
and  the  degree  of  fusion  is 

It  can  be  evaluated  once  A(q^  )  is  known.  The  results  for 
both  the  deuterium-tritium  mixture  and  the  tritium  gas  are 
shown  in  figure  12.  For  radiative  solution  is 

asymptotic  to  the  non-radiative  solution.  Away  from  the 
non-radiative  asymptotic  limit,  equation  (4.21)  depends  on 
the  specific  reaction  chemistry  through  A  which  is  a  function 
of  the  reaction  and  bremsstrahlung  rates. 

Certain  features  of  the  radiative  phase  space  solutions. 
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e.g.,  the  existence  of  extremal  points,  can  be  ascertained 
without  solving  the  radiative  phase  space  integral  equation. 
The  existence  of  extremal  points  on  the  shock  curve 

can  be  determined  by  differentiating  equation  (4.21)  with 
respect  to  ck  .  The  resulting  expression  is 


1  dA  _  128  /UJ  “  —  \  ^ 

HQ  doc  9  '  8  '  doC  ’ 


(4.22) 


Since 

JL_  M  =  R 

HQ  doc  p  ^ 

d<*.' 

solving  equation  (4.22)  for  T",  we  see  that 

d<c 


(4.23) 


^  _  _9_. 

d.c  ““  128 


£ 


(4.24) 


It  follows  that  uJ  has  an  extremal  value  when  the  bremsstrah- 

lung  power  loss  just  equals  the  thermonuclear  power  gain, 

(R  =1). 

P 

The  nature  of  the  extremal  point  of  equation  (4.21) 

depends  on  the  sign  of  the  second  derivative  evaluated  at 

the  extremum.  An  expression  for  — ^  is  obtained  by  dif- 

d<<2 

ferentiating  equation  (4.24)  with  respect  to  oC  .  The  result¬ 
ing  expression  is 


dJ; 

9 

'  (1-R  ) 

dc 

+  1 

dR 

d<2" 

128  ! 

d>< 

-1 

d->c 

(4.25) 


radiation  is  taken  into  account  in  the  structure  equations. 


the  speed  of  the  gas  increases  until  the  radiative  losses 

just  equal  the  gains  due  to  reactions.  The  speed  then 

begins  to  decrease  as  the  temperature  of  the  gas  falls  and 

the  density  of  the  gas  increases. 

If  the  state  equation  depends  on  oC  ,  as  it  does  for 

the  pure  tritium  case,  then  ~  eO  does  not  imply  that  ~  =  0. 

d*  d" 

In  fact,  it  follows  from  equation  (4.17)  that  when  Ai*-  =  0, 
d'^  —  TT 

then  -^= - ^jwhich  is  negative.  The  first  term  in  equa- 

d«  2(2+^) 

tion  (4.27)  is  K  .  Therefore,  if  is  negative, 

dR 

then  _ £  is  positive  and  the  extremum  is  a  maximum.  It 

dB 

can  be  shown  from  equation  (4.20)  that  is  negative  if 


2r  d(<y^'» 


(4.29) 


The  reaction  probability  for  T-T  reactions  has  the  form. 


<t\/ >=  const  X  (exp^-^^^j+ 


const  \  , 


(see  section  3.7).  Therefore  the  requirement  on  ~  ,  given 

in  equation  (4.29),  is  satisfied  and  the  extremal  point  on 

the  tritium  u  -oL  curve  is  a  maximum. 

The  existence  of  extrema  on  the  radiative  t-uC  curves 

can  also  be  determined  analytically.  For  the  deuterium-tritium 
gL  1 

case,  tT  =  0  implies  either  u;  =  —  or  “  =0.  When  the  entire 
’  d«<  2  d* 

system  of  equations  is  integrated,  we  see  that  the  dimension¬ 


less  speed  never  reaches  the  value  1/2.  Therefore  if  ^  =  0, 


then 


d** 


=  0  and  the  extrenal  values  of  temperature  and 


speed  occur  for  the  same  value  of  ^  ,  Since  l-2t0  is  posi¬ 
tive,  a  speed  maximum  corresponds  to  a  temperature  maximum. 

The  temperature  and  speed  extremal  values  do  not  occur 

at  the  seune  value  of  o(.  for  the  trititun  gas.  It  follows 

from  equation  (4.17)  that  vdien  ^  =  0,  then  ~  =  ”7^"  ■  .  >  0. 

^  doC  *  dd  2(l-2u)) 

Since  the  slope  of  theu-*  curve  is  positive  when  =0, 

da 

and  since  u:  goes  through  only  one  maximtun  value,  the 
temperature  maximum  occurs  for  a  smaller  value  of  rC  than 
the  speed  maximum.  In  fact,  the  physical  condition  for  the 
existence  of  a  temperature  maximum  is  not  that  the  bremsstrah- 
lung  power  equal  the  thermonuclear  power  as  it  is  for  the 
deuterium-tritium  case.  There  is  no  such  simple  criterion 
for  a  gas  whose  equation  of  state  depends  on  . 

The  coordinates  of  the  maximiun  points  on  the  radiative 
'X  -ci.  and  -  at  phase  space  curves  can  be  foimd  only  if 
the  phase  space  integral  equation  is  solved.  The  maximum 
on  the  deuterium-tritium  t  -  at  curve  occurs  at  the  point 
('t,  -X  )  =  (0.1019,  0.38)  •  The  maximum  on  the  U'  -at  curve 

occurs  for  (u,  at)  =  (0.285,  0.38)  .  A  comparison  between 
the  non-radiative  and  radiative  phase  space  solutions  shows 
that  the  effect  of  the  radiative  energy  loss  is  to  lower  the 
maximum  temperature  inside  the  wave  and  to  shift  its  occurr¬ 
ence  from  X.  =  0.89  to  X=  0.38  .  The  speed,  which 
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monotonically  increases  with  increasing  oC  when  radiation 
losses  are  neglected,  rises  to  a  maximum  value  and  then 
falls  when  radiation  losses  are  included  in  the  analysis. 

The  coordinates  of  the  maximum  point  on  the  ^  -  oC 
curve  for  tritium  are  ( 'ItT  >  ^  )  “  (0.0974,  0.37),  correspond¬ 
ing  to  a  dimensionless  speed  of  u)  =0.308.  Those  of  the 
maximum  point  on  the  lo  -  oC  are  (u)  ,  dC  )  =  (0.322,  0.56). 
When  radiation  is  included  in  the  equations,  the  temperature 
maximum  occurs  when  37%  of  the  material  has  fused  compared 
to  the  non-radiative  value  of  72%. 
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CHAPTER  V 

SHOCK  PROFILES  IN  REAL  SPACE 

5.1  Solution  of  the  Shock  Structure  Equation  for  a  Reacting 
and  Radiating  Gas 

The  profiles  of  the  temperature,  speed,  density,  pres¬ 
sure  and  degree  of  fusion  in  the  reaction  zone  can  be  found 
by  solving  the  full  system  of  shock  structure  equations.  The 
phase  space  solutions  are  used  to  eliminate  all  of  the  flow 
variables  except  d.  from  the  continuity  equation.  The 
resulting  non-linear  differential  equation  is 


where 

and 


doi 

dx  “  ^ 


(^)  , 


g  (o^ )  = 


^A  H 


<^«> 

10  2 


(5.1) 


^  =  0.6  ,  =  1, 


1  ,  1, 


♦t  = 


_ 1_ 

l+3<^ 


for  the  deuterium-tritium  mixture, 

for  pure  tritium,  and 

for  the  hydrogen-tritium  mixture. 


No  reactions  occur  in  front  of  the  wave.  This  physical 
fact  is  incorporated  into  the  differential  equation  by  assum¬ 
ing  that  the  reaction  probability  cuts  off  at  the  initial 
temperature  of  the  gas,  i.e.,  <lv'>(Tj^)=  0.  No  reactions 
occur  behind  the  wave  because  the  fuel  has  been  consumed,  i.e,o<.  = 
Therefore  g(oc)  approaches  0  as  x  approaches  +  »»  and  -  . 
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oC(x)  varies  directly  with  the  initial  number  density 
and  inversely  with  the  shock  sjjeed  since  AH  ~  .  This 

dependence  may  be  eliminated  if  the  continuity  equation  is 
made  dimensionless.  An  appropriate  scale  length  with  which 
to  make  the  continuity  equation  dimensionless  is  the  reaction 
mean  free  path  at  the  beginning  of  the  reaction  zone,  X  2* 

The  dimensionless  spatial  variable  becomes 


<5. 2) 


where  • 

The  dimensionless  continuity  equation  is: 


ai  -  Sah\2  (W  (5.3) 

ittp  W  ^ 

The  constant  AH  proportional  to  |n^/n2||<''Jy'Uj^. 

Since  n^^  and  n^  are  linearly  related  by  the  jump  conditions 

across  the  gas  dynamic  shock,  equation  (5.3)  is  independent 

of  the  initial  number  density  of  the  gas.  It  follows  from 

equations  (2.30)  and  (2.40)  that  is  proportional  to 

1/2 

T^  which  is  proportional  to  the  shock  speed.  The  different 

«• 

tial  equation  thereby  becomes  independent  of  also. 

In  most  shock  structure  problems,  one  must  integrate 
the  differential  equation  from  ^  =  -ooto  ^  =  +  0^ 
because  the  boundary  condition  is  known  only  at  ^  ^ 


tz  ^  CC 


1  m 


} 


!  i 


>  • 


r*'’  • 

yj 


>-rV% 
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In  this  problem,  the  boundary  condition  is  known  for  a  finite 
value  of  ^  which  for  convenience  we  assume  to  ^  =  0.  The 
assumption  of  the  von  Neumann  model  of  thermonuclear  shocks 
requires  that  oC  »o  throughout  the  ?as  dynamic  shock.  Al¬ 
though  ^  =0  ahead  of  the  gas  dynamic  shock,  “  ^0  behind  it 
d^  a? 

since  <'rv>(T2)  is  non-zero.  The  differential  equation  for 
^  (oC  )  in  the  reaction  zone  can  thus  be  reduced  to  the 
quadrature. 


(5.4) 


Equation  (5.4)  can  be  integrated  analytically  for 

o(.  2 

“<■<!.  UJ  and  <TV>  are  constant  over  the  range  of 
integration  and  can  be  taken  outside  the  integral.  Therefore, 

2  vt 


^  (  at  )  ^ 


"*T^2 


^AH>. 


(5.5) 


Performing  the  integration,  we  see  that 

1 


^  (X)  =  I!It 


2 

^2 


^AH  >2 


1  - 


- 1 


(5.6) 


This  expression  can  be  simplified  if  1/(1-  "Ji  )  is  expanded 
to  first  order  in  ^  .  The  resulting  linear  relation  be¬ 


tween  ^  and  <  is 
^  (X  )  ■ 


Sah.X2<‘*''4 


(5.7) 
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The  relation  between  ^  and  ck  ,  obtained  by  numerically 
integrating  equation  (5.4),  establishes  that  the  linear 

J 

relation  between  <5>  and  oC  is  valid  for  —  4  0.001. 

Equation  (5.4)  cazuiot  be  integrated  analytically  when 
2 

>  0.001.  uJ  and  remain  inside  the  integral.  Their 

functional  dependence  on  is  complicated.  A  closed  form 
expression  for  <<r\)>  (T: )  exists  only  for  T-T  reactions  at 
Mev  temperatures.  In  all  other  cases  must  be  evaluated 

numerically  (see  sections  3.2  and  3.7).  A  change  of  vari¬ 
ables  in  from  't  to  oC  results  in  an  integral  that 

cannot  be  evaluated  analytically.  Therefore  the  integration 
has  been  performed  namerically  with  an  Euler  forward  integra¬ 
tion  scheme. 

5.2  Solution  of  the  Shock  Structure  Ecmations  for  a  Radiating 
Gas 

The  shock  structure  equations  for  a  reacting  gas,  and 
for  a  reacting  and  radiating  gas,  were  presented  in  Chapter  IV. 
If  the  gas  radiates  but  does  not  react,  equations  (4.6)  and 
(4.14)  reduce  to 

2 

16  +  (  ^  +  2D)T  +  =  0,  (5.8) 

where  D  =  D*  +  4A, 

X 

and  ^  ^  ^  ^  • 
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i 
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4 

I 
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The  relation  between  the  dimensionless  speed,  r  >  and  A, 
obtained  from  the  energy  equation,  is 

i 

U)  »  (^2  CH  -  2A  -  10T)^  (5.9) 

An  iterative  method,  (similar  to  that  used  to  solve 
the  phase  space  equations  for  a  reacting  and  radiating 
gas),  can  be  used  to  solve  the  structure  equations  for  a 
radiating  gas ,  The  values  of  't  and  u.'  behind  the  gas 
dynamic  shock  (given  by  the  jump  conditions)  are  used 
for  the  first  iteration.  Since  these  values  of  and  oJ 
are  constant  for  all  x,  Aj^(x  )  (which  can  be  evaluated 
analytically),  varies  linearly  with  x.  Aj^(x)  is  substituted 
into  equation  (5.8)  to  obtain  '^^^(x).  Once  A^^  (x)  and  '^^^(x) 
are  known,  (x)  can  be  evaluated  with  equation  (5.9).  Given 
'Tj^(x)  and  UJj^(x),  A2(x)  can  be  computed  nvunerically .  The 
process  iG  repeated  until  '^^(x),  (x)  and  A(x)  remain 

the  same  to  some  required  accuracy  from  one  iteration  to 
the  next.  The  scheme  converges  to  three  decimal  places  in 
a  few  iterations  for  the  cases  we  have  considered. 

One  feature  of  the  no  reaction  limiting  case  of  the 
structure  equations  is  that  the  time-independent  flow  equa¬ 
tions  are  not  satisfied  for  all  values  of  x.  As  the  gas 
radiates,  its  temperature  and  speed  decrease.  7.*-  follows 
from  equation  (5.9),  that  for  values  of  HT  which  are  A4.  0,1, 


the  raaximiim  value  that  A  can  have  is  CH  =  0,5.  Therefore 
the  solution  of  these  equations  is  meaningful  only  for 
those  values  of  x  for  which  A  A.  CH.  The  cut  off  distance 
occurs  about  one  reaction  mean  free  path  behind  the  gas 
dynamic  discontinuity. 

5,3  Profiles  of  the  Flow  Variables  for  Shocks  Propagating 
with  the  Chapman -Jouguet  Speed  in  Deuterium-Tritium 
and  Tritium  Gases 

The  profiles  of  the  flow  variables  behind  the  gas 
dynamic  shock  have  been  computed  for  thermonuclear  shocks 
propagating  with  the  Chapman -Jouguet  speed.  The  profiles 
of  dimensionless  temperature,  speed,  density,  and  pressure 
are  shown  in  figures  16  to  19  for  the  deuterium-tritium 
mixture  and  figures  21  to  24  for  the  tritium  gas.  Three 
cases  are  displayed  in  each  figure:  (1)  pure  reaction, 

(2)  reaction  an'l  radiation,  (3)  pure  radiation.  The  pure 
reaction  curves  are  computed  for  values  of  o's.  up  to 
“C  =  0.99.  The  Euler  forward  integration  procedure  is 
inadequate  for  values  of  the  spatial  variable  corresponding 
to  rX.  greater  than  0.99.  However,  the  behavior  of  the 
curves  as  they  approach  the  values  given  by  the  jiimp  condi¬ 
tions  across  a  thermonuclear  shock  is  not  of  much  interest. 


The  reaction  and  radiation  curves  are  computed  for  values  of 
up  to  (X,  =0.80  for  the  deuterium-tritium  mixture,  and 
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c^  *  0.91  for  the  tritium  gas.  The  flow  variables  are  chang¬ 
ing  very  rapidly  with  space  for  larger  values  of  oC  ,  so  a 
higher  order  integration  scheme  is  needed.  However,  the 
behavior  of  the  variable  for  ^  >  0.80  is  of  no  physical 
interest.  The  pure  radiative  curves  are  computed  for 
values  of  x  such  that  the  radiation  energy  loss  pareuneter, 

A,  is  A  CH. 

The  reaction  and  radiation  shock  profiles  combine  the 
features  of  the  pure  reaction  and  pure  radiation  solutions. 

As  ^  ^  0,  all  three  curves  approach  the  values  given  by 

the  jximp  conditions  behind  the  gas  dyncunic  shock.  For 
^  ^  1,  the  reaction  and  radiation  solution  coincides 

with  the  pure  reaction  curve.  For  >  1,  it  follows  the 
general  shape  of  the  pure  radiation  curve.  However,  the 
two  do  not  coincide  because  the  composition  of  the  gas 
changes  as  the  reactions  occur,  and  the  pure  radiative  case 
is  valid  only  for  a  gas  consisting  of  reactants. 

The  pure  reaction  and  pure  radiation  shock  profiles  put 
upper  and  lower  bounds  on  the  values  of  the  flow  variables 
at  each  point  .  For  the  pure  reaction  case,  the  dimension¬ 
less  speed,  VA.  ,  rises  roonotonically  from  ^  =  0  to  ^  . 

The  dimensionless  density,  ^  »  ®nd  pressure,  t'  *»  1-  , 


decrease  monotonically  over  the  same  range.  The  temperature 


rises  to  a  peak  and  then  decreases  gradually  to  its  post¬ 
detonation  steady  state  value.  The  temperature  peak  occurs 
at  =13.9  for  the  deuterium-tritium  mixture  and  at 

=2.6  for  the  tritium  gas.  For  the  pure  radiation  case, 
all  of  the  variables  are  monotonic  functions  of  ^  .  The 
temperature  and  speed  of  the  gas  decrease  from  =  0  to 
^  =  oo  •  there  is  no  energy  source  (e.g.,  reactions)  to 

accelerate  and  heat  the  gas.  The  gas  density  and  pressure 
increase  monotonically  with  increasing  ^  . 

When  the  coupled  effects  of  reactions  and  radiation  are 
included  in  the  equations,  u.)  ,  4*  ,  and  IT  are  no  longer 
monotonic  functions  of  ^  .  At  first  the  speed  of  the  gas 
increases,  as  in  the  pure  reaction  case.  It  continues  to 
increase  until  the  power  input  due  to  reactions  is  equal  to 
the  power  loss  due  to  bremsstrahlung,  (R^  =  1) .  Then  the 
speed  begins  to  decrease,  as  in  the  limiting  pure  radiative 
case.  The  curves  of  dimensionless  density  and  pressure  also 
undergo  a  shift  in  slope  for  =  1;  the  pressure  and  density 
decrease  with  increasing  ^  for  R^  <  1  and  increase  with 
increasing  ^  for  R^  >  1.  The  turning  point  of  the  u)  ,  4  , 
and  "'T  curves  occurs  at  ^  =0.356  for  the  deuterium- 

tritium  case  and  at  =  0.947  for  the  tritium  case.  The 
temperature  increases  with  increasing  ^  until  R^  =  1  for 
the  deuterium-tritium  mixture  and  R  =  0.425  (  ^  =  0.386) 
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TEMPERATURE  PROFILE  BEHIND  GAS  DYNAMIC  SHOCK 
IN  DEUTERIUM  -  TRITIUM  MIXTURE  (hi) 
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FIG.  21  TEMPERATURE  PROFILE  BEHIND  GAS  DYNAMIC  SHOCK  IN  TRITIUM 


FIG.  22  SPEED  PROFILE  BEHIND  GAS  DYNAMIC  SHOCK  IN  TRITIUM 


BEHIND  GAS  DYNAMIC  SHOCK  IN  TRITIUM 


PRESSURE  PROFILE  BEHIND  GAS  DYNAMIC  SHOCK 
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25  DEGREE  OF  FUSION  BEHIND  GAS  DYNAMIC  SHOCK  IN  TRITIUM 
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for  the  tritium  gas.  It  then  falls  as  radiative  effects 
become  dcmiinant. 

The  variation  of  ,  the  degree  of  fusion,  with  space 
behind  the  wave  front  is  shown  in  figures  20  and  25  for  the 
two  gas  mixtures.  oC  increases  monotonically  with  increas¬ 
ing  ^  since  thermonuclear  reactions  are  irreversible. 

The  reaction  and  radiation  curve  coincides  with  the  pure 
reaction  curve  for  small  values  of  oC.  As  increases, 
oC  C  ^  )  increases  faster  for  the  reaction  and  radiation  case 
than  for  the  pure  reaction  case.  The  following  argument 
can  be  used  to  explain  this  result. 

The  ratio  of  the  reaction  rate  for  the  radiating  gas 
(R)  and  the  non -radiating  gas  (NR)  is 


NR 


fR 

2 

f  <1  >  (T  ) 

'nr  'nr' 


The  density  for  a  given  ^  for  a  radiating  gas  is  greater 
than  (or  equal  to)  thac  for  a  non-radiating  gas  (see  figures 
18  and  23).  The  temperature  for  a  given  ^  for  a  radiating 

I 

gas  is  less  than  (or  equal  to)  that  for  a  non-radiating  gas 
(see  figures  16  and  21) .  The  reaction  probability  for  T-T 
reactions  decreases  with  decreasing  temperature.  The  reac¬ 
tion  rate  is  proportional  to  the  reaction  probability  and 
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to  the  density  squared.  Since 

for  each  value  of  ,  the  reaction  rate  is  larger  (or  the 
Scuae)  vdien  radiation  losses  are  included  in  the  equations 
than  when  they  are  neglected  for  the  tritium  gas.  The 
reaction  probability  for  T-D  reactions,  (unlike  that  for 
T-T  reactions),  increases  with  decreasing  temperature  for 
kT>70  kev,  i.e,, 

<^''>  (V 

2  2 

Therefore,  even  if  =  p  , 

R  NR 

the  reaction  rate  is  greater  for  the  radiating  gas  than  it  is 

for  the  non-radiating  gas  for  the  deuterium-tritium  mixture. 

The  general  properties  of  the  shock  profiles  can  be 

disoissed  with  the  variardes  in  dimensionless  form,  but  in 

order  to  compute  the  thickness  of  the  reaction  zone,  the 

actual  values  of  the  physical  variables  must  be  kzu>wn.  For 

an  initial  temperature  of  5  x  10  ^K,  and  an  initial  pressure 

of  1.33  X  10^  dynes/an^,  (n^  10^^  (cm"^) ) ,  the  reaction  mean 

0 

free  path  behind  the  gas  dynamic  shock  is  ^  10  cm.  One 
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} 
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possible  definition  of  the  width  of  the  reaction  zone  is 
the  distance  over  which  80%  of  the  material  fuses,  A 
It  is  expressed  in  the  units  of  X  in  Table  5. 


Gas 

(kT) 2 (Mev) 

X2(cm) 

^.80 

pure  reaction  & 

reaction  radiation 

deuterivim- 

2.93 

6.01  X  10® 

5.43  1.54 

tritium 

tritiun. 

1.88 

1.59  X  10® 

4.83  2.96 

TABLE  5.  Reaction  Zone  Thickness 


The  distance  behind  the  gas  dynamic  shock  over  which  80% 
of  the  reactions  occur  is  the  order  of  several  earth  radii. 

The  thickness  of  the  reaction  zone  varies  inversely  with 
the  initial  pressure  of  the  gas.  Therefore,  it  will  be 
reduced  if  the  pressure  is  increased.  Reaction  zone  thick¬ 
nesses  for  smaller  degrees  of  fusion  can  be  computed  also. 
There  is  no  need  for  such  a  high  degree  of  fusion  in  order 
to  obtain  useful  power  from  thermonuclear  reactions. 

5.4  Profiles  of  the  Flow  Variables  for  Shocks  Propagating 

with  Chapinan-Jouquet  Speeds  in  Hydrogen-Tritium  Mixtures 

It  was  shown  in  section  2.3  that  the  minimum  shock  speed 

satisfying  the  Rankine*^ugoniot  conditions  for  a  non-radiative 
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thermonuclear  shock  is  the  Chapman- Jouguet  speed.  For 
shocks  travelling  with  the  C-J  speed,  the  fusion  reactions 
go  to  completion  inside  the  wave  and  a  steac^  state  is 

7 

established  behind  it.  Shock  speeds  of  the  order  of  10  cnv/sec 

are  presently  obtainable  in  electromagnetically  driven  shock 

0 

tubes.  Shock  speeds  of  10  cm/sec  have  been  reported.  There¬ 
fore,  a  study  of  thermonuclear  shocks  which  travel  with 
speeds  lower  than  the  C-J  speed  for  tritium,  or  a  deuterium- 
tritium  mixture,  is  of  interest  for  shock  tube  experiments. 

The  C-J  speed  for  a  hydrogen-tritium  mixture  is  a 
function  of  the  relative  concentration  of  the  two  species . 

Since  the  C-J  speed  for  the  mixture  is  proportional  to 

i 

,  (see  equations  (2.29)  and  (3.72)),  the  C-J  speed 

decreases  as  more  hydrogen  is  added  to  the  mixture.  The 

0 

C-J  speed  is  reduced  to  10  cm/sec  for  =  643,  Although 
the  power  yield  would  be  small  for  such  mj.xtures,  there 
would  be  a  measurable  neutron  flux  in  the  shock  tube.  Neu¬ 
tron  fluxes  can  be  measured  with  appropriate  counters.  Such 
measurements  could  be  used  to  determine  whether  the  plasma 
temperature  had  its  predicted  value  behind  the  gas  dynamic 
shock . 

The  flux  of  neutrons  per  unit  area  per  unit  time  is 
defined  as 


U 


(5.10) 


■ .  % 
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The  neutrons  are  assumed  to  be  In  equilibrium  with  the  other 
species  and  to  travel  with  the  local  fluid  speed.  If 
<<<■<.  — ,  it  follows  from  equations  (3,92)  and  (5.7) 


that 


**  <°‘y>2 


1 

(1+3  6 ) ^ 


X  . 


(5.11) 


The  distance j  x,  behind  the  wave  in  the  shock  reference 
frame  is  related  to  the  time,  t,  elapsing  after  the  wave 
has  passed  a  given  point  hy 


X  =  t  .  (5.12) 

Therefore  the  neutron  flux  also  increases  linearly  with  time 
behind  the  wave  for  •  The  values  of  the  neutron 

flux  at  2j  4,  6j  8  and  and  10  meters  behind  the  gas  dynamic 
shock  have  been  computed  for  hydrogen-tritivim  mixtures  with 

g 

(z  >*33.3  and  ^  »  643.  The  C~J  shock  speed  is  4.34  x  10 
cm/sec  for  mixtures  with  ^  =  33.3.  A  10  meter  gas  sample 
behind  the  gas  dynamic  shock  has  a  temperature  of  18.67  kev. 

g 

The  shock  speed  for  mixtures  with  6  =  643  is  1  x  10  cm/sec. 
The  temperature  behind  the  gas  dynamic  shock  is  0.981  kev. 
The  neutron  flux  as  a  function  of  space  and  time  is  given  in 
Table  6.  The  degree  of  fusion  corresjx^nding  to  each  value 


of  the  flux  is  also  included. 


V, 


7,  - 
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disadvantage  of  increasing  the  bremsstrahlung  rate  without 
increasing  the  reaction  rate.  Another  method  of  obtaining 
thermonuclear  shocks  might  be  to  send  shocks  with  speeds 
less  than  the  Chapman-Jouguet  speed  into  a  deuterium-tri¬ 
tium  mixture.  The  approximations  that  we  make  in  order  to 
solve  for  the  structure  of  shocks  travelling  with  the  C-J 
speed  become  better  for  shocks  travelling  at  lower  speeds 
(see  section  2.2). 

The  analysis  of  thermonuclear  shocks  propagating  with 

speeds  below  the  Chapman-Jouguet  speed  follows  very  simply 

from  the  analysis  of  the  C-J  case.  In  order  to  keep  the  jump 

conditions  across  the  gas  dynamic  shock  in  the  same  form 

as  those  for  the  C-J  case,  we  define  a  constant  Q^,  such 

that  the  C-J  relation  between  the  shock  speed  and  Q  is 

m 

preserved  for  any  value  of  the  shock  speed,  i.e., 

u?  =  2(  y  ^  -  1)  .  (5.13) 

1  m 

Then,  §  .  (5.14) 

2 

The  shock  structure  equations  for  deuterium-tritium  mixture, 

Qm 

presented  in  section  3.4,  are  valid  only  for  oC<  ix  =  —  , 

m  Q 

where  Q  is  the  value  of  Q  when  the  shock  speed  is  equal  to 


the  C-J  speed.  The  time  independent  equations  are  not 
satisfied  for  larger  values  of 
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Temperature  profiles  in  the  reaction  zone  have  been 

0  0 

computed  for  shock  speeds  of  1.10  x  10  cnv/sec,  3.48  x  10 

9 

cnv/sec  and  1,10  x  10  cm/sec  (which  correspond  to  o(  ^  =  0.001, 
0.01,  and  0.1,  respectively).  The  results  are  displayed 
in  figure  26.  The  profiles  of  the  relative  degree  of  fusion. 


'V 


,  in  the  reaction  zone  are  shown  in  figure  27.  A  graph 


of  the  ratio  of  the  bremsstrahlung  power  to  the  thermonuclear 
power,  R^,  for  the  reaction  zone  is  shown  in  figures  28  for 
the  three  cases  and  figure  29  for  shocks  travelling  with  the 
C-J  speed  in  a  deuterium-tritium  mixture  and  a  trititun  gas. 

Radiative  effects  are  not  important  for  shocks  travel- 
ling  with  speeds  between  3.48  x  10  and  1.10  x  10  cnv/sec. 

-3 

R^  remains  at  the  roughly  constant  value  of  6  x  10  for 

0 

shocks  travelling  with  a  speed  of  3.48  x  10  cn/sec  and 

-2  9 

2.5  X  10  for  shocks  propagating  at  1.10  x  10  ar/sec.  The 

temperature  profiles  for  the  two  cases  look  like  the  non- 

radiative  temperature  profile  for  shocks  travelling  with  the 

C-J  speed;  the  temperature  rises  to  a  peak  value  and  then 

decreases  as  approaches  .  The  details  of  the  shape  of 

m 

each  curve  depend  on  the  temperatures  occurring  in  the 
reaction  zone  (which  depend  on  the  shock  speed) .  The  reaction 
probability  does  not  have  the  same  functional  dependence  on 
temperature  for  all  temperatures.  The  temperature  and  relative 
degree  of  fusion  rise  fastest  for  temperatures  close  to  the 


TEMPERATURE  BEHIND  GAS  DYNAMIC  SHOCK  IN  DEUTERIUM  -  TRITIUM  MIXTURE 


FIG.  28  RATIO  OF  RADIATIVE  POWER  Lv,.3  TO  THERMONUCLEAR  POWER  BEHIND  GAS 
DYNAMIC  SHOCK  IN  DEUTERIUM -TRITIUM  MIXTURE  (1:1)  FOR  DIFFERENT  SHOCK  SPEEDS 


RATIO  OF  RADIATIVE  POWER  LOSS  TO  THERMONUCLEAR  POWER  BEHIND  GAS  DYNAMIC 

SHOCK 
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naxixaum  value  o£  the  reaction  probability  (kT^70  kev) . 

Radiative  effects  becorae  important  as  the  shock  speed 

0 

is  reduced  to  1.10  x  10  cn/sec.  The  power  ratio  in  the 
reaction  zone  is  8  x  10  The  larger  power  ratio 
occurs  because  the  reaction  probability  decreases  rapidly 
as  kT  falls  from  40  kev  to  4  kev.  The  power  ratio  decreases 
with  dist2mce  behind  the  gas  dynamic  shock.  Even  though 
fuel  is  being  replaced  by  inert  reaction  products  which 
radiate  but  do  not  react,  the  reaction  probabilility  is 
rising  very  fast  for  the  temperatures  occurring  in  the 
reaction  zone  (3-4  kev)  so  that  the  reaction  rate  is  increas¬ 
ing  faster  than  the  bremsstrahlung  rate.  The  temperature 
profile  does  not  look  like  the  other  profiles  shown  in 
figure  26  because  radiative  effects  are  significant.  The 
temperature  increases  with  distance  behind  the  gas  dynamic 
shock  but  does  not  reach  a  peak  and  then  decrease  because 
the  condition  for  a  temperature  maximum,  Rp"l>  is  never 
satisfied. 

One  definition  of  the  thickness  of  the  reaction  zone 
is  the  distance  from  the  beginning  of  the  reaction  zone 
to  the  point  where  a  fixed  percentage  of  the  gas  has  reacted. 
The  reaction  zone  tliicknesses  for  shocks  travelling  with 
the  three  different  shock  speeds  given  above  have  been  computed 
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for  four  different  percentages  of  reacted  gas.  The  results 
are  presented  in  Table  7. 


#  cm  . 
u,  (  —  ) 

1  sec 

Partial  Thickness  of  Reaction  Zone  (cm) 

0.01% 

0.1% 

1% 

10% 

0.001 

1.10  X  10® 

8.5  X  10^ 

7.9  X  10® 

0.01 

3.48  X  10® 

7.0  X  10^ 

7.3  X  10® 

1.4  X  10® 

0.1 

1.10  X  10® 

2.7  X  10® 

2.7  X  10^ 

2.9  X  10® 

7.9  X  10® 

1.0 

3.48  X  10® 

2.8  X  10^ 

2.8  X  10® 

2.8  X  10® 

3.3  X  10^ 

TABLE  7.  Thickness  of  Reaction  Zone  for  a  Given  Degree 
of  Reaction  as  a  Function  of  Shock  Speed 

The  thickness  of  the  reaction  zone  depends  on  the  tempera¬ 
tures  and  densities  occurring  in  the  reaction  zone  and  the 
shock  speed.  rises  to  a  maximum  value  for  kT  70  kev 

and  then  falls  off  with  increasing  temperature.  One  would 
think  that  the  reaction  zone  thickness  would  be  minimized 
for  shock  speeds  such  that  =  70  kev.  However,  the 
thickness  depends  on  shock  speed  and  the  square  of  the  density 
as  well  as  on  reaction  probability.  The  thicknesses  of  the 
reaction  zone  (for  small  values  of  <a  )  can  be  compared  for 
two  shocks  travelling  with  different  shock  speeds  in  the 


.1 


I 


-  r 


following  way.  it  follows  from  the  continuity  equation  that 


\  /a 


ado(\  Ax. 

/  ^  •  Since  the  initial 

a  ^  a  2 

density  of  the  gas  is  the  same  for  both  cases,  is 

2  Av 

~  .  The  ratio  of  the  mass  flow  constant,  ,  is 

/  NT  ^ 

equal  to  /  1  .  Therefore, 


If  T^  is  >  T^  ,  then 


. 


is  <  1 


r.b  v'ir 


(For  example,  if  T^  =  70  kev  8.9  x  10  cm^/sec) 

and  T2  =  29.3  kev  (4<rv^'5:  6.6  x  lO”^^  cmVsec  ),  then 


^  1.1  .) 


One  can  determine  the  shock  speed 


that  gives  the  thinnest  reaction  zone  from  the  graph  of 
reaction  probability  vs.  temperature  for  T-D  reactions.  The 


minimum  thickness  for  a  given  degree  of  reaction  occurs  for 


/V^ 


40  kev  which  corresponds  to  a  shock  speed  of  4.08  x  10 
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The  gas  dyneunic  shock  was  assuined  to  be  r;aich  thlmier 

than  the  reaction  zone  for  shocks  travelling  at  the  C-J 

speed.  The  validity  of  this  assximption  for  shocks  tra^^• 

elling  with  speeds  below  the  C-J  speed  will  now  be  examined. 

The  thickness  of  the  gas  dynamic  shock  is  assumed  to  be  one 

collision  mean  free  path  evaluated  at  the  end  of  the  gas 

dyncunic  shock.  It  is  found  to  be  168  m  for  shocks  propagating 
0 

at  3.48  X  10  oiv/sec.  (The  collision  mean  free  path  varies 
roughly  as  the  temperature  squared.)  Therefore,  although 
it  takes  only  73  m  for  0.1%  of  the  gas  to  react  for  a  shock 
speed  of  3.48  x  10  cnv/sec,  168  m  or  more  are  required  for 
the  gas  temperature  to  build  up  to  its  final  value  behind 
the  gas  dynamic  shock. 

The  behavior  of  the  flow  variables  over  a  distance  of 

10  m  behind  the  gas  dyneunic  shock  is  relevant  to  shock  tube 

experiments.  The  neutron  flux  as  a  function  of  space  (or 

time,  if  viewed  in  the  laboratory),  is  shown  in  Table  8 

0  0 

for  shock  speeds  of  1.10  x  10  av'sec  and  3.48  x  10  cn/sec. 

The  corresponding  values  of  ot.  are  also  noted.  Much 
higher  neutron  fluxes  are  obtainable  from  deuterium-tritium 
mixtures  than  from  hydrogen-tritium  mixtures  for  shocks 
of  the  same  speed. 

For  thermonuclear  shocks  to  be  obtainable  in  shock 
tubes,  the  thickness  of  the  gas  dynamic  shock  must  be  minimized 


4"  '•^t*^-*  ^  ♦‘•r-eT*3i«^'s!Sfe 
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However,  the  shock  speed  must  be  high  enough  so  that  reac¬ 
tion  rates  are  appreciable.  The  gas  dynamic  shock  is 

1.68  m  thick  for  a  thermonuclear  shock  propagating  at 

8  17  *2  “I 

1.10x10  cm/sec.  Neutron  fluxes  of  10  (cra'^-sec”'^) 

are  achieved  in  the  reaction  zone.  About 


TABLE  8.  Neutron  Flux  Behind  Gas  Dynamic  Shock  in  a 
Deuterium-Tritium  Mixture  as  a  Function  of 
Shock  Speed 
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0.0001%  of  the  gas  will  have  reacted  8  m  behind  the  gas 
dynamic  shock.  The  specific  thermonuclear  power  is 


TN  D2  T2  ^  '2  R 


.-16.2 . -“18.  _ _  _ _  .-“13. watts 

(2  X  10  )  (1.66  X  10  )  (17.6  x  1.602  x  10  )— rr— 

cm** 


2  kw/cm  . 


If  the  shock  tube  is  10  m  long  and  has  a  radius  of  10  cm, 
then  the  total  power  output  is  cs  5  x  10^  kw.  (Fission  power 

reactors  have  power  outputs  of  the  same  order  of  magnitude.) 

A  shock  tube  is  a  pulsed  device  so  the  power  output  is 
not  continuous.  A  estimate  can  be  made  of  the  total  energy 
released  in  reactions  with  each  pulse.  We  assume  that  a 
given  gas  element  will  react  from  the  time  the  shock  wave 
initiates  the  reactions  until  the  rarefaction  wave  following 
the  shock  wave  quenches  them.  If  the  shock  wave  has  travelled 
for  a  time,  t,  the  distance  between  the  shock  front  and  the 
rarefaction  wave  (in  the  limit  of  very  large  initial  Mach 


number)  is  (Gross  ) 


Ul 

Xit)  -  T 


(5.15) 


The  energy  released  during  the  operation  of  the  shock  tube 


may  be  estimated  by 


Ht)  dt 


)  "I 


2>  111  !o£ 


(5.16) 
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The  operating  time,  is  equal  to  the  length  of  the 

apparatus  divided  by  the  shock  speed.  If  the  shock  speed  is 
8 

^  10  aa/aec,  then  the  operating  time  for  a  10  m  long 
shock  tiibe  is  c:  10  yusec.  Evaluating  equation  (5.16), 
(noting  that  2kw/cm^  and  r=10  cm) ,  we  see  that  the 

energy  released  is  c;  750  joules.  This  particular  figure  is 
of  course  only  a  rough  estimate.  However,  it  indicates 
that  an  energy  output  of  the  order  of  hundreds  of  joules 


is  feasible 
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constant  that  depends  on  A,  B  and  C 
enthalpy 

Planck's  constant 

mass  emission  coefficient 

Boltzmann  constant 

mass  absorption  coefficient 

constant  in  shock  structure  differential  equation 

de  Broglie  wavelength 

Debye  shielding  distance 

photon  absorption  mean  free  path 

reaction  mean  free  path 

reaction  mean  free  path  behind  gas  dynamic  shock 

alpha  particle  mass 

deuteron  mass 

electron  mass 

neutron  mass 

proton  mass 

triton  mass 

effective  triton  mass 

alpha  particle  number  density 

deuteron  number  density 

electron  number  density 

neutron  number  density 

proton  number  density 
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n 

Q 


Q 


eff 


q 


R 


r 

<r 
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^Gov 
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triton  number  density 

dimensionless  speed 

power  due  to  bremsstrahlung 

power  due  to  thermonuclear  reactions 

pressure 

dimensionless  pressure 
dimensionless  density 

energy  released  in  reactions  p/er  gram  of  initial 
gas  mixture 

effective  energy  released  per  gram  of  initial 
gas  mixture 

energy  released  per  reaction 
radiative  energy  flux 
thermonuclear  energy  per  unit  mass 
universal  gas  constant 

ratio  of  bremsstrahlung  power  to  thermonuclear 
power 

density 

cross  section 

Gamov  cross  section 

Govorov  c<.oSS  section 

reaction  probability 

temperature 

dimensionless  temperature 


self  collision  time 


Tee 


Tr 


u 

W 

w. 

J 


self  collision  time  for  electrons 
self  collision  time  for  neutrons 

equilibration  time 
reaction  time 
gas  speed 

shock  speed  (initial  gas  speed) 
mean  molecular  weight 
production  rate  of  ith  specie 
dimensionless  spatial  variable 


z 


atomic  number 
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